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Chapter  1 
INTRODUCTION 


Electromagnetic  transmission  through  apertures  in  a  conducting 
screen  of  finite  thickness  has  been  studied  extensively  by  many  in¬ 
vestigators.  The  effect  of  screen  thickness  on  the  transmission  is 
of  interest  in  many  cases.  The  simplest  model  of  the  problem  is  a 
two-dimensional  one  and  is  the  focus  of  most  investigations.  Several 
solutions  have  been  developed  for  the  problem  of  electromagnetic  pene¬ 
tration  through  an  infinitely  long  slit  in  a  thick  screen  [1]— [10] . 

Both  TM  and  TE  cases  have  been  studied,  and  a  variety  of  techniques 
are  used.  However,  most  investigators  treat  only  the  problem  of  a  slit 
with  rectangular  cross  section.  Morita  [3]  and  Auckland  [10]  developed 
solutions  for  two-dimensional  slits  with  arbitrary  cross  sections. 

Uslenghi  [23]  also  considered  the  problem  of  a  two-dimensional  gasket 
in  a  screen  under  several  simplifying  assumptions. 

A  three-dimensional  problem  is  of  greater  complexity.  Host  of  the 
related  work  has  been  done  for  the  quasistatic  case.  Akhiezer  [11]  extended 
Bethe's  [12]  theory  to  the  case  of  a  circular  aperture  in  a  thick  screen. 
Garb  examined  the  problem  of  a  narrow  rectangular  slot  in  a  screen  of 
finite  but  small  thickness  [13] ;  and  later  studied  the  polarizability  of 
small  openings  in  a  thick  screen  [14],  McDonald  [15]  considered  the 
polarizabilities  of  small  circular  and  rectangular  apertures  in  a  thick 
screen.  The  methods  of  solutions  usually  rely  on  the  assumption  that 
the  aperture  is  electrically  small,  and  also  that  the  waveguide  modes 
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are  known  for  the  aperture  region  inside  the  screen. 

This  report  considers  the  more  general  problem  for  which  an 
aperture  in  a  thick  screen  is  rotationally  symmetric  about  an  axis 
perpendicular  to  the  screen.  The  method  of  solution,  in  general, 
does  not  depend  on  the  size  of  the  aperture.  The  formulation  is  not 
limited  to  the  situation  where  the  aperture  region  in  the  screen  is 
a  familiar  waveguide  section  for  which  the  characteristic  modes  are 
known.  However,  this  case  is  also  considered  and  investigated  in 
detail. 

In  Chapter  2,  we  develop  the  basic  formulation  of  the  problem. 

The  equivalence  principle  [16,  Sec.  3-5]  is  used  to  separate  the 
problem  Into  different  regions  of  interests  by  means  of  unknown  boundary 
currents.  Field  operators  and  operator  equations  are  then  established 
for  the  currents.  A  modal  formulation  is  used  to  treat  the  special  case 
when  the  aperture  region  is  a  waveguide  region  (in  our  case  a  cylin¬ 
drical  or  coaxial  region) .  In  the  more  general  case  a  nonmodal  formula¬ 
tion  Is  used.  Chapter  3  contains  a  moment  method  [17]  approach  to  the  solu¬ 
tion  of  the  simultaneous  operator  equations.  Both  modal  and  nonmodal 
formulations  are  given  and  the  measurement  matrix  is  discussed.  A  low 
frequency  analysis  is  presented  in  Chapter  4  for  a  filled  narrow  annular 
slot.  Resonant  behavior  of  the  power  transmission  is  observed,  and  the 
transmission  coefficient  and  the  electric  polarizability  are  discussed. 
Numerical  results  are  presented  in  Chapter  5  for  several  examples. 

Emphasis  is  put  on  the  electric  and  magnetic  currents  obtained.  Modal 
and  nonmodal  results  are  compared  when  both  formulations  apply.  Also, 


numerical  results  and  analytical  predictions  for  low  frequency  power 
transmission  are  compared.  A  final  discussion  is  presented  together 
with  some  recommendations  in  Chapter  6. 
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Chapter  2 

PROBLEM  FORMULATION 


2.1.  Problem  Specification 

The  problem  to  be  considered  is  shown  in  Fig.  1,  which  shows 
a  conducting  screen  of  thickness  d  situated  between  the  surfaces  z  =  0 
and  z  *  d.  The  regions  z  £  d  (denoted  region  c)  and  z  £  0  (denoted 
region  a)  are  coupled  through  an  aperture  region  (denoted  region  b) 
in  the  screen.  The  incident  fields  and  H^,  generated  in  region  a, 
penetrate  through  the  aperture  region  and  radiate  into  region  c. 
and  H  are  the  fields  produced  by  the  impressed  sources  and  iL  in 
unbounded  space  filled  with  the  same  material  as  in  region  a.  The 
aperture  region,  region  b,  is  bounded  by  the  following  three  surfaces 

the  interface  between  regions  a  and  b  which  lies  in 
the  plane  z  =  0. 

82^  the  interface  between  regions  b  and  c  which  lies  in  the 
plane  z  =  d. 

S^:  the  interface  between  region  b  and  the  conducting  screen. 

S^,  and  are  all  assumed  to  be  rotationally  symmetric  about  the 

z-axis.  The  three  regions  are  assumed  to  contain  media  characterized 

by  their  permittivities  and  permeabilities,  (Cg,  Ua)  >  (G^»  U^)  and 

(t  ,  u  ) .  For  simplicity,  all  three  regions  are  assumed  lossless, 
c  c 

A  cross  section  of  this  situation  is  shown  in  Fig.  2(a). 


For  the  special  case  where  region  b  is  a  cylindrical  waveguide 
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Fig.  1.  Transmission  through  a  rotationally  symmetric  aperture  in 


a  thick  conducting  screen. 
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region,  a  conductor  in  the  center  of  this  region  can  support  a  TEM 
mode  and  thus  increase  the  power  transmission.  Therefore,  we  gen¬ 
eralize  our  problem  of  investigation  to  one  for  which  the  boundary 
surface  can  be  either  simply-connected  (as  shown  in  Fig.  2(a))  or 
doubly-connected  (as  shown  in  Fig.  2(b)).  The  boundary  surfaces 
and  are  then  either  circular  (when  is  simply-connected)  or 
annualar  (when  is  doubly-connected).  Since  S^,  S2  and  are  all 
symmetric  about  the  z-axis,  their  union,  S  =  S^uS^uS^,  is  also  sym¬ 
metric  about  the  axis.  In  other  words,  S^,  S^,  and  their  union 
are  all  surfaces  of  revolution.  A  generating  curve  F,  as  well  as  its 
associated  coordinate  t,  and  unit  tangent  vector  j:,  can  be  defined  for 
S.  Note  that  F  is  open  with  its  end  points  on  the  z-axis  when  is 
simply-connected  and  closes  upon  itself  when  S^  is  doubly-connected. 

F  is  the  union  of  the  generating  curves  of  S^,  and  S^,  denoted  F^,  F^ 
and  F^.  A  unit  normal  vector  is  defined  for  every  point  on  S  as 

n  =  *  1  (2-1) • 

A 

Here  $  is  the  conventional  unit  vector  associated  with  the  azimuth  angle 
<f>.  The  arrangement  described  above  is  illustrated  in  Figs.  3(a)  and  3(b) 
for  the  two  possible  situations.  Furthermore,  each  boundary  surface  can 
be  approached  from  either  of  the  two  regions  that  the  surface  separates. 

We  denote  the  side  of  each  surface  from  which  the  normal  vector  n  points 
by  the  superscript  "+",  and  the  opposite  side  by  the  superscript  For 

example,  according  to  the  convention  we  use  in  this  work,  as  shown  in 


Fig.  3(a).  Generating  curves  of  S^,  S2  and  and  the  unit 
vectors,  simply  connected. 


r 


Fig.  3(b).  Generating  curves  of  S^,  S2  and  and  the  unit 
vectors,  S„  doubly  ..’nnected. 
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Figs.  3(a)  and  3(b),  S+  denotes  the  surface  immediately  next  to  just 
outside  of  region  b. 

2.2.  Equivalent  Problems 

In  this  section,  the  equivalent  principle  [16,  Sec.  3-5]  is  used 
to  separate  the  original  problem  into  three  equivalent  situations  which 
exist  in  the  three  regions  a,  b,  and  c.  For  region  a,  a  magnetic  cur¬ 
rent  sheet  M  is  placed  just  inside  the  conducting  plane  at  z  =  0,  with 
covered  with  conductor  also.  The  equivalent  magnetic  current  is 
defined  as 


M.  =  £  x  E  (2-2) 

where  E  is  the  electric  field  over  z  =  0  in  the  original  problem  and 
is  unknown  on  S^.  It  has  zero  tangential  (to  the  z  =  0  plane)  com¬ 
ponent  elsewhere  in  the  z  =  0  plane.  This  magnetic  current  M^,  and  the  im¬ 
pressed  sources  and  M^,  radiate  in  the  presence  of  the  complete 
conducting  screen  over  z  =  0  to  give  the  correct  fields  E  and  H  in 
region  a.  This  situation  is  shown  in  Fig.  4(a).  Furthermore,  image 
sources  can  be  used  to  account  for  the  effect  of  the  infinite  conduct¬ 
ing  plane  as  far  as  filled  in  region  a  are  concerned.  Therefore,  we 
have  the  equivalent  situation,  shown  in  Fig.  4(b),  established  for 
region  a.  Here  J^,  M^,  and  M^,  together  with  their  images,  radiate  into 
unbounded  space  filled  with  a  medium  characterized  by  (e  ,  y  )  to  create 
the  correct  fields  in  region  a.  A  similar  equivalent  situation  can  be 
developed  for  region  c  as  shown  in  Fig.  5.  Here  the  equivalent  magnetic 
current  defined  by 


H 

—a 


<e 


V 


conducting  plane 


% 


Fig.  4(a).  Equivalent  situation  for  region  a. 


z=0 


;i 


E  ,  H 
—a  —a 


v; 


(images  of 


Fig.  4(b).  Equivalent  situation  for  region  a,  using 
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M2  =  K  x  2  ,  (2-3) 

together  with  its  mirror  image  with  respect  to  the  plane  z  =  d,  exist 
on  S2  and  radiate  into  unbounded  space  filled  with  the  medium  char¬ 
acterized  by  (£c>  1*c)*  These  currents  give  the  correct  fields  and 
in  region  c.  Note  that  in  the  two  equivalent  situations  for  regions 
a  and  c,  electromagnetic  fields  from  the  equivalent  magnetic  currents 
can  be  found  using  the  f  ield  operators  for  unbounded  space. 

For  region  b,  the  entire  closed  boundary  S  is  replaced  by  a 
perfect  electric  conductor,  and  equivalent  magnetic  current  sheets  -M 
and  -M2  are  placed  just  inside  and  S2>  respectively.  This  situation 

is  shown  in  Fig.  6,  where  -M.  and  -M^  radiate  in  the  presence  of  the 

closed  conductor  to  give  the  correct  fields  E,  and  H,  in  region  b.  Note 

— b  b 

that  the  use  of  -M^  and  -M^  in  this  region  ensures  the  continuity  of  tan¬ 
gential  electric  field  across  and  S2-  In  the  special  case  where 
region  b  is  a  cylindrical  or  coaxial  waveguide  region,  fields  can  be  ex¬ 
pressed  in  terms  of  and  M2  by  means  of  the  waveguide  modes.  However 
in  a  more  general  case,  it  is  not  possible  to  obtain  the  fields  in 
region  b  directly  from  and  M2  alone.  Therefore,  for  the  general  case, 
we  replace  the  conducting  surface  S  in  Fig.  6  by  the  electric  current, 
denoted  -J,  induced  on  S  by  -M^  and  •  This  equivalent  situation  is 
shown  in  Fig.  7,  where  -M^ ,  -M2>  and  -J  radiate  into  unbounded  space 

filled  with  a  medium  characterized  by  (c,  ,  |i,  ).  They  give  zero  fields 

b  b 

just  outside  S  and  the  correct  fields  E,  and  H,  in  region  b.  Note  again 

— b  — b 

that  we  have  established  an  equivalent  situation  in  which  only  the  field 
operators  for  unbounded  space  are  needed.  To  summarize,  we  have 
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developed  equivalent  situations  shown  in  Figs.  4(b)  and  5  for  regions  a 
and  c  respectively.  The  equivalent  situation  shown  in  Fig.  6  is  intended 
for  the  special  case  where  region  b  represents  a  waveguide  region,  while 
the  equivalent  situation  shown  in  Fig.  7  is  suited  to  the  general  case. 

In  the  special  case  where  region  b  is  a  waveguide  region,  if  waveguide 
modes  are  used  for  the  field  operators  needed  in  the  equivalent  situ¬ 
ation  in  Fig.  6,  we  call  the  formulation  a  modal  formulation.  If  a  gen¬ 
eral  treatment  corresponding  to  the  equivalent  situation  in  Fig.  7  is  used, 
we  call  the  formulation  a  nonmodal  formulation. 


2.3.  Field  Operators 


From  the  discussion  in  the  previous  section,  it  is  evident  that  two 
types  of  field  operators  are  needed  for  our  formulation.  These  are  the 
potential  integrals  which  give  the  electric  field  or  magnetic  field  in  an 
unbounded  medium  due  to  an  electric  or  magnetic  surface  current  distribu¬ 
tion,  and  the  modal  representation  which  gives  us  the  fields  inside  a 
waveguide  region  when  the  tangential  electric  field  is  specified  at  the 
two  end  surfaces.  In  the  first  case,  although  both  the  electric  and 

magnetic  types  of  currents  and  fields  are  considered,  because  of  duality 

6  h 

[16,  Sec.  3-2],  only  the  following  two  basic  operators,  and  are 
needed : 


^(J) 


-  iK 


J(r')e 


-jk  R 
J  n 


4ttR 


da'  +  -rf-  V  , 
jk  J 
J  a  1 


(V’  •  J(r ’ ) )e 


-jk  R 

J  a 


4ttR 


da'  (2-4) 


and 


o[(— )  *  V  x 


J(r')e 


-jk  R 


4ttR 


da' 


(2-5) 
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The  integrations  in  both  (2-4)  and  (2-5)  are  over  the  surface  A  where 
the  surface  distribution  J  resides,  and 


■/ 


p*  +  p'^  -  2pp'  cos  (<£'  -  4>)  +  (z 


’.'Y 


(2-6) 


n. 

k  =  /u>  y  e 

a  J  a  a 


(2-7) 


If  r_  approaches  A  from  one  side  of  the  surface,  and  if  the  surface  is 
assumed  to  be  smooth  there,  then  we  can  write  [21] 


Lh(J)  -  \  J (r)  xfi  -fix 
a  —  l  —  -  o  — o 


n^x  j(r')  x  (r-r')  (1  +  jk  R)e 

3 

4ttR 


-jk  R 
J  a 


da' 


+  n  [n  •  L  (J)] 
— O  — o  — 


(2-8) 


where  n  is  the  unit  normal  vector  of  the  surface  at  r  which  points 
— o  - 

toward  the  side  of  the  surface  from  which  r  approaches  A.  Note  that 

J. 

the  component  of  L^(J)  that  is  normal  to  A  is  not  written  explicitly 
because  it  is  not  used  in  our  operator  equations.  It  follows  that 
for  a  surface  electric  current  distribution  J  and  a  surface  magnetic 
current  M  in  an  unbounded  medium  characterized  by  (t^,  y^) ,  the  electro¬ 
magnetic  fields  E(J,  M)  and  H(J,  M)  are  given  by 


E(J,  M)  =  Le(n  J)  -  Lh(M) 

- —  (x  ot—  a  ~ 

H(J,  M)  =  ~  [L®(M)  +  Lh(naJ)] 


where 


(2-9) 

(2-10) 

(2-11) 


Note  that  the  linearity  of  the  operators  is  used  in  (2-9)  and  (2-10). 
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To  discuss  the  modal  type  of  operators  we  consider  the  case  in  which 
we  have  a  waveguide  region  b  between  z  =  0  and  z  =  d  formed  by  two  con¬ 
centric  cylindrical  conductors,  located  at  p  =  and  p  =*  Rout-  The 
boundary  conditions  are,  from  (2-2),  (2-3)  and  the  equivalent  situation 
shown  in  Fig.  6: 


and 


z  x  E,  =  M  at  z  =  0 

—  — b  1 


(2-12) 


z  X  E,  =  -  M  at  z  =  d 

—  — D  ""Z 


(2-13) 


Also  we  have  the  condition  that  the  tangential  electric  field  is  zero  on  the 

waveguide  walls.  Note  that  when  R  =  0  the  problem  reduces  to  one  in 

which  is  simply  connected  and  the  cylinder  in  the  center  is  missing. 

To  establish  the  operators  in  this  case  we  first  assume  that  the  E^  and 

H,  are  generated  by  two  vector  potentials,  A  =  A  z  and  F  =  F  z  [16,  Sec. 
b  z  z 

3-12].  A  is  called  the  magnetic  vector  potential  and  £  is  called  the 
electric  vector  potential.  The  electric  and  magnetic  fields  they  generate 
are: 


4  " 


and 


jkn^A  +  —  V(V  -  A)  -  V  x  F 
b  fr- 


-ikv 

F  +  — -  V(\7*  F)  +  V  X  A 

nb  "  Jkbnb 


(2-14) 


(2-15) 


Az  and  satisfy  the  Helmholtz  equations 

72A  +  k2A  =0  (2-16) 

z  b  z 


2  2 
V  F  +  RTF  = 
z  b  z 


0 


(2-17) 


and  the  fields  they  generate  satisfy  the  boundary  conditions  on  the 

waveguide  walls.  can  be  expanded  as  a  linear  combination  of  the 

TEM  modes,  and  A__w,  and  a  set  of  TM  modes,  {A  ,  A  n=0,±l,±2,.. 

l  km  l  km  nm  nm 

m=l, 2, 3 , . . . } .  F^  can  be  expanded  as  a  linear  combination  of  a  set  of  TE 

modes,  {F+  ,  F  I  n=0, ±1, +2, . . . ,  m=l,2,3, . . . } .  The  waveguide  modes  are 
nm  nm 

defined  in  the  following: 

Z'(k,  ,z)  (defined  only  for  R,  ^  0)  (2- 

b  in 


A1  -  zYk  _,z) 


nm  nm 


nm 


(2- 


(2- 


where 


f  N  (x  )J  tr-2-  x  )  -  J  (x  )N  (— ^ —  x  ) 
l  nnmnR  nm  nnmnR  nm 
'  out 


out 


if  R.  i  0 
in 


1 

n 


J  (x  ) 

2(1  -  4 — — ) 


4m(p) 

nm 


J„(Yx  ) 
n  nm 


(2-21) 


L 


/2  J  (-1—  x  ) 
n  R  nm 

out 


x  |  j  .  (x  Tf 

nm  n+1  nm 


if  R.  =  0 
in 


N’(x’  )J  (— ~  x’  )  -  J’(x’  )N  (— £—  x'  ) 
n  nm  n  R  4  nm  n  nm  n  R  ^  nm 
out  out 


2  J,2(x'  ) 

2 [1  ~  -  (1  -  ■7n-r,)  %  — 1 
x*  Y  x ' Z  J’Z(Yx'  ) 


if  R,  /  0 
in 


nm  n  ntn 


Z  ~  (k,  z)  =  e~]K  if  k  i  0 


Z  (0,z)  =  z 


(2-23a) 

(2-23b) 


Z  (0,z)  =  1 


(2-23c) 


and 


Y  = 


in 


(2-24) 


out 


K  -  ^2 

out 


if 


V- 


>  x 
out  —  nm 


nm 


J  -  A 

V  out 


if  k  R  "  x 
n  out  nm 


(2-25) 


x 


!<  ~  2 

out 


if  k,  R  >  x' 
d  out  —  n 


nm 


k'  =  ' 
nm 


x 


.  /  ,  nm  n2  ,2 

H 

Y  out 


if  k,  R  <  x' 
b  out  nm 


J  ,  N  ,  J',  N'  are  the  nth  order  Bessel  functions  of  the  first  and  second 
n  n  n  n 


kind,  and  their  respective  derivatives,  x  and  x'  are  the  mth  lowest 

nm  nm 


positive  real  roots  of  the  following  equations: 


N  (x  )J  (Yx  )  -  J  (x  )N  (Yx  )  =  0 
nnmn  nm  nnmn  nm 


if  R.  #0  (2-26) 

in 


N'(x’  )J'(Yx'  )  -  J'(x’  ) N '  (Y x '  )  =  0 
n  nm  n  nm  n  nm  n  nm 


J  (x  )  -  0 
n  nm 


if  R,  =0  (2-27) 

in 


J'(x’  )  =  0 
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Note  that  the  TEM  modes  can  only  exist  when  R.  ^  0.  Now  we  examine 

in 

the  electric  and  magnetic  fields  associated  with  each  of  the  modes. 


—  —  /s 

To  compute  the  fields,  we  substitute  and  A  z  in  (2-14)  and  (2-15) 


TEM± 


,  _m  t  , 
and  E  ,  and 
— nm 


for  A  to  obtain  their  corresponding  electric  fields  £' 

TEM±  m+  + 

their  magnetic  fields  H  and  H  .  We  substitute  F  z  into  (2-14) 

—  — nm  - 


and  (2-15)  for  £  to  obtain  the  corresponding  fields 
results  are  the  following: 


nnr~ 
,ex 


+  e+ 

E“  and  H  . 
nm  — nm 


The 


„TEM" 


TEM  +^kbZ 
5  e 


(2-28a) 


nm 


jn4>±ik  z  Y 
m  nm  nm 

e  e  — ; — — 

— nm  t  jwe. 


c-niB-A  ^  + 

'r  ’  k  z 
out  nm~ 


if  k  M 
nm 


(2-28b) 


m+  m 

E  «  e  e 
— nm  — nm 


Em 
— nm 


jn<l>  +  kj  A+  z 


nm— 


A  z 
nm— 


if  k  =0 
nm 


(2-28c) 


Ee±  -  ee  ejn4>  Z”  (k*  ,  z) 
nm  — nm  nm 


(2-28d) 


htem1=  etem± 

+  nb  -  - 


,m ' 


r.m: 


H  =  r  Y  z  *  E 
— nm  +  nm  —  — nm 


if  k  /  0 
nm 


(2-29a) 

(2-29b) 


„•+  .  .  j  li!  i  *  E"+ 

— nm  n,  — 

D 


Hm_  =  -  z  x  em  ein* 
— nm  —  — nm 


if  k  =0 
nm 


(2-29c) 


e+  .  e  +  1 

H  -  t  Y1  z  x  Ee  +  ( 

-nm  +  nm-  nm  jtop. 


~-)2  F  z 
R  „  nm— 
out 


if  k’  t  0 
nm 


(2-29d) 
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i 

•  i 
i 

v  i 


1 

Mr* 


,.e+  1  ~  ..  „e-f  .  „-f-  ^ 

H  - - : - : -  z  x  h  -  Joit,  F  z 

— nm  i  rj,  k,  z  —  —run  b  nm— 

b  b 


He_  «  F~  z 
— nm  nm— 


if  k'  «  0 
nm 


(2 


where 


TEM 


yin<Rout/Ri,) 


ni  m  /  v  *\  in  . m  /  \ 

e  *  \p  (p)p  +  (p)  (p 

— nm  rnm  —  p  nm  x 


e 

e 

--nm 


^  ^m(p>p  +  r  (p>£ 

P  nm  —  nm 


(2 


(2 


(2 


i  f  x  <  k  R 

nm  —  b  out 


nm 


1  _  (_-0HL_)2 

Vout 


/,  nm  ,2 

inb  /  (kTR~  -  1 

\/  b  out 


1  _  ( — EHH_y 

kKR  b 

b  out 


(2 


if  x  >  k.R 

nm  b  out 


Y' 

nm 


if  x'  <  kLR 

nm  —  b  out 


(2 


i  U^-)2  ~  1 


L 


V. 


out 


if  x'  >  k,  R 
nm  b  out 


The  set  of  sectors  {eTEM,  em  ,  e6  |m  =  1,2,3,...)  for  a  fixed  n  is  an 

—  — nm  — nm 

orthonormal  set  under  the  inner  product: 


:-29e) 

-30) 

-31) 

-32) 

-33) 

-34) 
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m 

«*5 


(a,  b)  =  a*  •  b  pdp 


C2-35) 


We  now  write  the  electric  field  E,  in  region  b  as  a  linear  combination 

b 


of  all  possible  modes: 


+  jV  -  "jkb\  TEM 


E,  =  (E,  •  z)z  +  (<i  e  +  ft  e  )  e 

b  — b  -  o  o  — 


+  l  7  [ 8+  Z+(k’  |Z>  +  8  Z  (k’  ,z))  ee  ejn4, 

^  L-  nm  nm  — nm 


n=-<A  m=l 


nm  nm 


'*>  a'  ,  ik  z  -jk  z  .  , 

c  r  .  +  nm  ,  -  nm  m  jncp 

+  )  /  («  e  +  u  e  )  e  e 

u  ,  nm  nm  — nm 

n=-°°  m=l 

k  /  0 
nm 


V  r  +  m  1  nth 

+  )  )  a  e  e 

L'  L.  nm  --nm 

n=-tr'  m=l 

k  =0 
nm 

From  (2-28),  (2-29)  and  (2-36),  we  can  write  for  H,  : 

~ b 

H  ^  1  ,  +  jRbZ  +  -  -jkb\  ,  TEM 

H,  =  (H,  •  z)z  +  —  (-rt  e  +  a  e  )  h 

— b  “b  “  —  r\  o  o  — 

D 

^  jk  z  -jk  z  .  . 

,  r  r  .,  /  +  nm  ,  nm  .  ,  m  jnd) 

+  >  )  Y  (- 1  e  +  a  e  )  h  eJ  Y 

u  ,  nm  nm  nm  -mm 

n=-rt'  m=l 


k  *0 
nm 


i  k  *  z  i\.  ** 

v  r  ,, ,  ,  .,+  nm  ,  n-  ,‘'nm‘\  ,e  jnd) 
+  )  )  Y  (-jf,  e  +  8  e  )  h  eJ  T 

'■  u  nm  nm  nm  — nm 


-jk'  z 


n=-”  m=l 

k'  *  0 
nm 


nm  nm 


+  8  e  )  h  e 

nm  — nm 


+  )  /  — (k  z  a  +  B  )  h  eJ  ^ 

-  n,  b  nm  nm  — nm 
n=—1r  m=l  b 

k  =0 

nm 


+  >  l  — V  e+  i'e  ejn^ 

L  L  -  -  j  rit  k,  nm  ^nm 
n=-fY  m=l  b  b 

kf  =0 
nm 


(2-36) 


1 


(2-37) 
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where  B  is  a  constant  and 
nm 


TEM 


z  x  e 


TEM 


hm 
— nm 


a  rn 

z  x  e  , 
—  — nm 


,e  '  „  e 

h  *  z  x  e 
-nm  —  — nm 


(2-38) 


From  (2-12),  (2-13),  (2-36)  and  the  orthonormal  relationships  among 


TEM  m  e  , 

e  ,  e  and  e  ,  we  obtain 
-  — nm  — nm 


%  +  =  (^TEM’ 


+  jv  +  -  ~ikbd 

a  e  +  a  e 
o  o 


(hTEM,  M°) 


+  .  ,,m  wn. 

a  +  <x  =  (h  ,  M, ) 
nm  nm  -nm  — ] 


,  jk  d  -]k  d 

. +  „  nm  nm 

a  e  +  ft  e 

nm  nm 


-  (h™  ,  M") 
—nm  —2 


for  k  ^  0 
nm 


+  ,,  m  n 

a  =  •  h  ,  M, 
nm  — nm  —1 


-  (h"m,  M") 


for  k  =0 
nm 


B+  +  0"  =  (he  ,  M?) 

nm  nm  — nm  — ] 


,  jk'd  - jk '  d 

nm  ,-  nm 

p  e  +  p  e 

nm  nm 


-  <*w  *S> 


for  k'  ^0 
nm 


3  -  (he  ,  m?) 

nm  — nm  —1 


B+  d  +  B 

nm  nm 


-  (^m’  M2> 


for  k'  =  0 
nm 


(2-39) 


where 


wn  .,n  a  l  wn  i  1  .  - 

^  ~  Mip& +  V- =  s 


or 


2  7r  2tt 

(M.  •  p)e-jnit>  d4>  +  £  (M1  •  ^)e"jn4>  d<t>] 

0 


^i  =  1  (Mif£  +  MV)C 

n=-«'  ' 


J 

0 

i  n<J> 


i  =  1,2 

i  =  1,2, 


(2-40a) 

(2-40b) 
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From  (2-39) ,  we  obtain 


.  ,.TEM  +jkbd  wo  ,  .,o. 

±  ±j(h  ,  e  M:  +  M  ) 

ao  2  sin  k,  d 

D 

+  ,  -  ,,  TEM  XA°\  /.TEM  U°  .  ,  j  „0. 

aQ  +  aQ  =  (h  »  M^),  (h  ,  +  cos  k^d  M^) 


if  k  d  ^  TT,2ir, . . . 
b 


(2-41a) 


k^d  =  tt,2tt,  . . . 


+ j  k  d 

■  .  /.in  nm  ..n  ,  ..n. 

±  j  (\m>  e  +  V 


2  sin  k  d 
nm 


if  k  d  ^  0,tt, 2rr, .  . . 
nm 


(2-41b) 


a+  +  a  =  (hm  ,  M"),  (hm  ,  m"  +  cos  k  d  m")  =  0  if  k  d=ir,2iT, . .  . 
nm  nm  — nm  —1  — nm  —1  nm  —2  nm 


+  ,,  m  wn.  ,,m  wn  ,  ,,n.  _ 

a  =  (h  Mi )  >  (h  _>  Mi  +  M,)  =  0 

nm  — nm  —1  nm  —1  —2 


if  k  d  =  0 
nm 


+j  k'  d 

.  .  e  nm  n  n. 

J(hnm>  e  — 1  +  — 2^ 

2  sin  k '  d 
nm 


B+  +  (T  =  (h®  ,  M"),  (h6  ,M"+cosk'  d  m") 
nm  nm  — nm  -1  — nm  — 1  nm  —2 

-(h6  ,  M?  +  M?)  p 

0  =  - 2"L_^ - —  ,  ff  =  (h  ,  M") 

nm  d  nm  nm  —1 


if  k'  d  ^  0,tt,2tt,  . . . 
nm 

(2-41c) 

if  k'  d  =  it  ,  2tt  , . . . 
nm 


if  k'  d  =  0 
nm 


We  substitute  (2-41)  into  (2-36)  and  obtain 


(Efe  •  £)z  +  e 


-esc  k^d  (h^^,sinkb(z-d)M°+sinkbzM2)  if  sink^d  ^  0 
(hTEM,  M^coskkZ  +  j  (ot^-aQ) sink^z  if  sink^d = 0 


}  }  esc  k  d(hm  ,  sink  (z-d)M^1  +  sink  zM^)  e™  e 

L  '■  nm  — nm  nm  —1  nm  — 2  — nm 

n=-m  m=l 

sink  d^O 
nm 


J  }  csck  d  (h  ,  sink  (z-d)M1+sink  zM-)e  e 
L  '■  nm  — nm  nm  —1  nm  —2  — nm 

n=-on  m=l 

sink'  d^O 
nm 
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oo  oo 


+  J  J  [  (h  ,  M  )cosk  z  +  j  (et  -a  )smk  z]e  e 
^  — nm  —1  nm  J  nm  nm  nm  nm 

n=-<»  m-l 


m  jn4> 


k  d=7T,27T,... 
nm 


00  oo 


+  I  I  [  Ov  ,M  )cosk  z  +  j(B  -8  )sink'  z] 
L  L.  1 nm  -1'  nm  J  nm  nm  nm 

n=-oo  m=i 


e  Jn( 4> 

e  e 
— nm 


kf  d=TT,  2tt,  .  .  . 
nm 


oo  oo 


,  V  r  .rn\  m  Jn(J' 

+  )  >  (h  ,  M, )e  e 

L  L ,  — nm  — 1  —nm 

n=-oo  m=i 


k  =0 
nm 


00  OO 


+  l  y  (hV.  (1  -  X  -  f  M?)e 

-  ~nm  cl  — 1  a  — z  — 

n=-oo  m=l 


— \  m"  _  Z  \  eJn(J> 

— nm 


kr  =0 
nm 


We  substitute  (2-41)  into  (2-37)  and  obtain 


— b  ’  Lb(~-1’  -a2> 


TEM 


:sck^d(  h  ,  coskhCz-d^^+cosk^zM^) 


J  ,  TEM 

(Hb  *  z)z  -  h 

^  i  (hTKM,Mt.1)sink,  z  +  C  cosk,z 

— 1  b  o  b 


7  y  j  Y  esck  d(h  ,cosk  (z-d)M1  +  cosk  zM_)h 
-  nm  nm  — nm  nm  —1  nm  —2  — 

n=-oo  m=l 


,n'hm 
— nm 


sink  d/0 
nm 


y  />  7  Y  csck  d(h  ,cosk  (z-d)M,  +  cosk  zM„)h 

*- ,  nm  nm  — nm  nm  —1  nm  —2  — i 

n=-oo  m=]_ 


n'he 
— nm 


sink'  d^O 
nm 


•O  00 

I  V 


l  i  Y  [ (h  ,H. ) sink 
*•,  nm  -nm  —1  i 

n=-«  m=l 

k  d=7T,27I,... 
nm 


+-  C  cosk  z]hm  e^n<^ 
nm  nm  — nm 


(2-42) 

if  sink^d  j*  0 
if  sink^d  =  0 
ejn<J) 

ejn(t> 
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e  jn<J> 


J  J  Y'  h  ,M  )sink  z  +  D  cosk  z]h  e 
^  ,  nm  — nm  — 1  nm  nm  nm  — nm 

n=-oo  m=l 


k'  d*n,2iT,.., 
run 


l  l  ^  .M?)k.z  +  B  Ih"1  eJn<J> 

L  ln,  — nm  -1  b  nm  — nm 

n**-00  m=l  b 


k  =0 
nm 


™  .  (he  ,  m”  +  M")  .  . 

V  y  x  '-nrc  -1  -2  he  eJn4> 

^  r>  k,  d  -nm 

n=-oo  m=l  b  d 

k'  =  0 


(2-43) 


where 


C  =  -  j  (a  -  a),C 
o  J  o  o  nm 


j(a  -  -X  ) ,  D  *=  -  j(3  -8  )  (2-44) 

nm  nm  nm  nm  nm 


w 

L,  (-M.  ,  -M  )  gives  the  magnetic  field  inside  the  waveguide  region  b  for 
D  ~1  ~~L 

the  equivalent  situation  in  Fig.  6.  Note  again  that  the  TEM  mode  can 
only  exist  when  /  0.  Also,  we  note  that  in  (2-43),  corresponding  to 
each  undetermined  constant  (not  determined  from  and  ^  but  can  be  de¬ 
termined  through  other  information)  C  ,  B  ,  C  or  D  ,  there  is  a  con- 

°  o  nm’  nm  nm 

dition  on  the  magnetic  currents  as  shown  in  (2-41).  This  correspondence 
is  summarized  in  Table  1.  The  z-romponents  of  the  fields  in  the  above 
analysis  are  not  investigated  in  detail  because  they  are  not  used  in 
the  operator  equations  which  we  introduce  in  the  next  section. 
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Situation 

Undetermined 

Coefficient 

Condition 

on  the  Current 

11 

2t r, . . . 

c 

o 

(hTEM, 

Mj)  = 

.  ,  ..TEH 

-  cosk^d  (h  ,M2) 

ii 

o 

B 

(hm  , 

M?)  - 

-  (hm  ,  m") 

nm 

nm 

— nm 

—1 

— nm  —2 

k  d  =  tt. 

2tt,  . .  . 

c 

(hm  , 

M?)  = 

-  cosk  d(hm  ,M"l 

nm 

nm 

— nm 

—1 

nm  — nm  2 

k'  d  =  tt. 

2ir.  — 

D 

(h6  , 

M?)  = 

-  cosk'  d(h6  ,M?) 

nm 

nm 

-nm 

—1 

nm  -nm  —2 

Table  1.  Summary  of  correspondence  between  undetermined  coefficients  in 
w 

L  and  conditions  on  the  magnetic  currents, 
b 


2.4.  Basic  Operator  Equations 

So  far  we  have  developed  the  equivalent  situations  for  the  three 
regions  a,  b  and  c.  In  each  region  the  electromagnetic  fields  can  be 
obtained  from  the  surface  equivalent  currents  through  the  operators  in¬ 
troduced  in  the  last  section.  For  region  a,  magnetic  current  2M^  on  S^, 
together  with  the  impressed  sources  and  and  their  images  J_^  and  M^, 

radiate  into  unbounded  medium  (f  ,  M  )  to  give  E  ,  H  ,  as  shown  in 

a  a  — a  —a 

Fig.  4(b).  Therefore,  we  have 


E 

-a 


n-Xii +  vi » -  l>i  + «;»  -  2La(v 


(2-  45' 


—a  ■  n"  +  «;>  +  '-a<Vl  +  Vi»  +  f  La%> 

a  a 


(2-46) 


where  the  field  point  is  anywhere  in  the  region  z  <  0.  For  region  c. 


0  V  .  1  Ifcl  *k'-  .  m  ■  — -  -1-*-  -  -  T 


the  situation  is  very  similar  except  that  no  impressed  sources  are  in 
this  region.  Therefore,  we  have 


For  region  b,  we  have  magnetic  currents  -M  on  S.  and  -M„  on  S„.  These 

— 1  1  — z  Z 

magnetic  currents,  together  with  the  electric  current  on  the  entire  S, 

radiate  into  unbounded  medium  (£,,  g,)  to  give  E,  ,  H,  in  region  b.  There- 

b  b  — b  — b 

fore,  we  have 


h  '  -  lb(V>  +  +  >{%> 


(2-49) 


\  '  -  £  lLb%>  +  +  <2-50) 

b 

If  we  substitute  (2-4)  and  (2-8)  into  (2-45)  and  (2-47),  we  find  that 
tangential  E  is  zero  everywhere  on  the  plane  z  =  0  except  in  S-.  .  In 

3  -L 

S^,  let  the  field  point  approach  from  z  <  0  (i.e.,  n^  =  -  £  in  (2-8)) 

and  obtain  £  x  =  M^.  Also  £  x  E^  is  zero  everywhere  on  z  =  d  except 

in  S2,  let  the  field  point  approach  from  z  >  d  (i.e.,  n^  =  £  in  (2-8)), 
and  obtain  £  x  E^  =  -  Mj.  These  boundary  conditions  are  expected  because 
of  the  way  the  equivalent  situations  were  set  up  for  regions  a  and  c.  For 
region  b,  however,  the  condition  that  fields  be  zero  outside  S  still  has 
to  be  enforced.  This  can  be  done  by  requiring  either  the  tangential  elec¬ 
tric  field  or  the  tangential  magnetic  field  be  zero  just  outside  S.  We 

choose  to  use  the  condition  on  the  electric  field. 
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If  (2-51)  is  satisfied,  we  immediately  have  from  (2-49)  that  z  *  E. 

—  — b 

is  zero  on  and 


z><E=M=zxE 
—  — b  — 1  —  — a 


zxe,  =  -M=zxe 
—  — b  —  2  —  — c 


(2-52) 


(2-53) 


Therefore,  if  (2-51)  is  satisfied,  our  tangential  electric  field  will 
be  zero  on  all  conductor  surfaces  and  continuous  through  the  two  aper¬ 
ture  faces  and  S2.  The  boundary  condition  yet  to  be  satisfied  is 
the  continuity  of  tangential  magnetic  field  through  and  S^.  To 

compute  tangential  H,  on  S.  and  S„,  we  note  that  when  (2-51)  is  satisfied 

b  1  L 

Ampere's  law  states 


nXnxH,  =  -nxJ  in  S  and  S 

—  —  ~D  JL  Z 


(2-54) 


To  compute  the  tangential  in  S ^ ,  we  note  that  the  contribution  from 

the  impressed  sources  and  their  images  is  twice  that  radiated  by  J^, 

into  unbounded  (e  ,  y  ).  Therefore 
a  a 


nxnxH  =  nxnx  [2H,  +  — -  L  (Mn  )  ] 
—  —  —a  —  —  —I  T)  a  —1 

a 

The  tangential  in  S2  is  simply,  from  (2-48), 


2  6 

n  x  n  x  H  =  n  x  n  x  —  l  (M„)  on  S„ 
-  -  -c  -  -  nc  c  -2  . 


on  S, 


(2-55) 


(2-56) 


Therefore,  the  continuity  of  tangential  magnetic  field  can  be  written  in 


the  following  form: 
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n  x  n  x  Le(M )  +  n  x  j  =-n  n  x  fi  x  h  on  S.  (2-57) 

3  J.  Z  Si  1  1 

„  „  e  n0  , 

n  x  n  x  Lc(M2)  +-^fixj=0  on  S2  (2-58) 

Equations  (2-57),  (2-58),  and  (2-51)  rewritten  below, 

n  x  n  x  [I^Q^)  +  l£(M2)  -  ^(n^J)]  =0  on  S+  (2-59) 

form  the  basic  operator  equations  for  the  general  problem,  where 

M^,  M2  and  .J  are  the  unknowns. 

When  region  b  represents  a  waveguide  region,  the  fields  in  region 
b  can  be  written  in  terms  of  and  M2  alone.  The  electric  current  J 
is  no  longer  needed  and  the  boundary  conditions  on  the  waveguide  walls 
are  built  into  the  waveguide  modes.  Therefore,  all  we  need  to  consider 
is  the  continuity  of  the  tangential  magnetic  field  across  and  S2: 

n 

n  x  n  x  Le(M,)  +  ~  n  x  a  x  l*(M,  ,  Mj  =  -  r)  n  x  fi  x  h  on  S,  (2-60) 
si  4  b  1  4  3  i  -L 

n 

n  x  fi  x  L®(M2)  +  -f  n  x  fi  x  M2)  =0  on  S2  (2-61) 

(2-60)  and  (2-61)  are  the  basic  operator  equations  for  the  special  case 
when  region  b  is  a  waveguide  region. 

2.5.  Fourier  Decomposition 

So  far  we  have  not  taken  advantage  of  the  rotational  symmetry  of 
the  geometry  assumed  in  the  problem.  In  this  section  this  property  is 
used  to  reduce  the  system  of  equations  to  a  smaller  set  of  equations. 

To  do  this,  we  first  define  the  Fourier  coefficients  of  the  unknowns 


as  follows: 
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Jn(t)  -  j"(t)£  +  jJ(t)$. 
2tt 


t 

2ir  J 


i 


2tr 


[J( t ,4>>  •  t]e  ^  d<{>  +  — 


[J(t,4>)  •i]e"Jn<t>  d4>  (2-62a) 


Mj(t)  =  Mjt(t)£  + 
2tt 


2tr 


2^  {  •  £]e-^n^  d<j)  +  —  [M^t,*)  .f]e“jn<}>  d4»  (2-62b) 

0  0 


M"(t)  =  M"t(t)£  +  M^Ct)! 


2tt 


2tt 


•  £3e_;intt  d4)  +  j  [M2 (t ,4»)  •ile~jnv  d<f>  (2-62c) 

6  0 

where  n  =  0,  ±1,  ±2,...  .  Since  the  field  point  in  the  operator  equations 
is  always  on  the  surface  S  and  only  the  tangential  component  is  used,  it  is 
convenient  to  define  field  quantities  in  the  same  fashion  as  for  the  sources 
in  (2-62): 


4 


2tt 


C(J)  " 


£  r2lT  $  211 

[L ®<J)  •  |]e~jn<})  d<»  +  f  [L®(J)  •  £]e"Jn<tl  d0 


0 


0 


2tt 


Lhn±(j)  =  — 

a  ^  2tt 


.  hn+ 


n£<j>-iie-i-*d*+4' 

0 


2ir 


[Lj(J)  •ile'jn4>  d<t> 


(2-63) 


(2-64) 


where  for  ,  the  field  point  is  on  the  side  of  the  surface  that  n  points 
away  from.  L^n  in  the  meanwhile,  is  used  for  field  point  on  the  other  side. 
Substituting  (2-4)  and  (2-8)  into  (2-63)  and  (2-64),  we  find 


Len(J)  =  Len(Jn) 

a  —  a  — 

=  [Len  (j”)  +  Len.  (j")]t  +  [L6?  (Jn)  +  L6”  (j”)]$ 
att  t  ott<t»  4>  -  cx4>t  t  v 


(2-65) 
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where 


Lh"  (f) 

a4>4> 


jk3 

a 

2v 


(z'-z)G3f  (t’)pMt' 


cos  v  =  £ 


z  ,  sin 


v  =  e 


•  p 


cos  v'=  £*•  z  ,  sin  v'=  £'•  p 


gn(t,f)  = 


-jk  R 

a 


cos  n<{> 


k  R 

a 


dcf) 


G^t.t')  « 


-jk  R 

T(1+JkaR)e  a  .2* 

- — x -  sin  -b  cos  nq>  d4> 

k  R 

t  a 


G2(t,t')  = 


r  a  +  jkR)e 

i  \  _  i  r< 


-jk  R 

a 


3  3 

.  k  R 
0  a 


cos  $  cos  ncf)  d<j> 


G^  (t , t ' )  = 


-jk  R 

1T  (1  +  jk  R)e  a 
J  a  _ 

3  3 

k  R 


sin  0  cos  n(J>  d<J> 


R  = 


/ 


(p  -  p*)  +  (z  -  zf)  +  4ppf  sin^ 


£ 

2 


( 2— 67h) 


(2-68a) 

(2-68b) 


(2-69) 


(2-70) 


(2-71) 


(2-72) 


(2-73) 


Note  that  in  (2-67),  f(t)  represents  the  source  function  and  all  integrals 

are  along  the  generating  curve  where  f(t)  is  not  zero.  For  the  operator 
w 

L,  ,  if  we  define 
b 


p 

2tt 


2tt 


tLb(-r  V 


£]e  d(J) 


2tt 


^>1 


m2) 


£]e  d<J> 


(2-74) 


it  follows  from  (2-43)  that 


csckj}d(hTE^,cosk^(z-d)M°+coskj32M2)hTE:M  if  n=0,  and  sink^d^O 


,j  r  .  TEM  o .  .  ,  .TEM 

q  l (h  .M^isink^z  +  C^cosk^zjh 


if  n=0,  Rin^0  and  sink^d=0 


if  n=^0  or  R.  =0 
In 


+  i  )  Y  esc  k  d(h  ,  cos  k  (z-d)M^  +  cos  k  zM_)h 

J  *•,  nm  nm  — nm  nm  —1  run  —  2  — 

m=l 


n'hm 
— nm 


sink  d^O 
am 


+  J  l 

m'l 


Y'  esc  k'  d(he  ,  cos  k'  (z-d)M^1  +  cos  k1  zM^)h6 
nm  nm  -nm  nm  —1  nm  — 2  -nm 


sink'  d^O 
nm 


+  i 


m=i 


Y  [(hm  ,  M?)  sin  k  z  +  C  cos  k  z]hm 
nm  —nm  -1  nm  nm  nm  — nm 


k  d-1!,!-, 

nm 


+  j  /. 

m=l 


Y'  [(he  ,  M?)  sin  k'  z  +  D  cos  k'  z]he 
nm  nm  —1  nm  nm  nm  — nm 


k'  d  =tt  ,  2  rr  . . 
nm 


+  J- 

\ 


[  (hm  ,  m")  k,  z  +  B  ]hm 
1  -1  nq  — nq 


-nq 


(h!_,  M"  +  m") 

_zm_rL  2 


kbd 


h 

-nq 


if  there  exists  k  =0 
nq 


if  there  exists  k'  =0 
nq 


(2-75 
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Note  that  all  x  and  x'  for  a  fixed  n  are  interlaced,  and  therefore 
nm  nm 

at  most  one  of  all  k  and  k'  can  be  zero.  Also  note  that  the  con- 

nm  nm 

ditions  on  the  currents  as  shown  in  Table  1  should  be  considered  as 
part  of  the  definition  of  the  operator. 

From  (2-64)  to  (2-67),  we  can  rewrite  (2-57)  to  (2-59)  as: 

n 

-j  en/wnN  .  cL^  n  _,n  „ 

-L  (M, )  +  n  x  J  =  n  H.  on  T.  (2-76) 

a  —1  l  —  —  a— l  1 

n 

-Len(M”)  +  n  x  Jn  =  0  on  F  (2-77) 

c  —  2  l  —  —  l 


Thn+,..n.  ,  thn+.  n.  Ten.  ,n.  _ 

-Lb  (m^)  +  (m^)  -  (n^  J  )  -  o  on  r 


(2-78) 


From  (2-74)  and  (2-75),  we  can  rewrite  (2-60)  and  (2-61)  as: 

n 


en.  n.  a  wn.  n  n.  _  ,.n 

"La  (-V  -  T  La  (^1>  —2  =  na^i 


on  r. 


(2-79) 


-L“nQg)  M°)  =  0 


on  r„ 


(2-80) 


For  both  sets  of  equations  above,  we  have  n  =  0,  ±1,  ±2,.. 
defined  as 


.  and  is 


h!  p  +  h"  : 

-1  lp- 


V 


2tt 


2tt  I  'Si 


p)e  d(|»  + 


A 

277 


277 


(Hi  *  i)e  jn<i> 


(2-81) 


Our  problem  has  now  been  reduced  to  solving  (2-76)  to  (2-78)  for  the 
general  case  and  solving  (2-79)  and  (2-80)  for  a  special  case.  Although 
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(  an  infinite  number  of  modes  exist,  higher  order  modes  are  not  important 

unless  frequency  is  high.  Also,  for  axial  plane  wave  incidence,  H,  is 

— i 

constant  in  and  therefore,  from  (2-81),  only  the  n  +  1  modes  are 
excited. 
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Chapter  3 


NUMERICAL  SOLUTION 


3.1.  Generating  Curve,  Basis  Functions  and  Symmetric  Product 

The  geometry  of  our  problem  can  be  fully  specified  by  the 
generating  curve  T.  This  curve,  in  general,  can  be  fully  described 
by  a  pair  of  parametric  functions 

p  =  h1(t),  z  =  h2(t)  (4>  fixed)  (3-1) 

In  our  numerical  solution,  a  finite  number  of  points,  (p  ,  z^) 

(p  ,  z2),...,  are  specified  on  the  curve.  T  is  approximated  by 
connecting  successive  points  with  straight  line  segments.  t^  de¬ 
notes  the  t-coordinate  of  the  ith  point  (p^,  z^) .  Therefore,  we  have 


p  =  hl(t)  =  p.  +  : >  (pi+1  -  pa) 

l 


2  =  h2(t)  =  2i  +  (zi+i  ~  zi) 


where 


At  =  t .  -  t . 

i  l+l  l 


(3-2) 


(3-3) 


(3-4) 


For  r„,  we  have  t.  <  t  <  t„  .  For  F, ,  we  have  t„  <  t  <  t„  .  For  T  , 

2  1  —  —  N^  1  Np  —  —  3 

if  S^  is  doubly  connected  (if  there  is  a  center  conductor  in  region  b) , 

we  have  t„  <  t  <  t  and  t„  <  t  <  t  ,  . . .  If  S.  is  singly  connected, 
N,  —  —  N„  N_  —  —  N.+l  3 

1  2  3  4 

we  only  have  t  £  t  £  t  for  F  .  This  arrangement  is  shown  in  Fig. 

N1  N2  J 
8(a)  and  Fig.  8(b)  for  the  two  cases. 
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We  now  define  a  set  of  basis  functions,  (u  . ,  u, . I  i  •  1,2,... N  ; 

— ti  -=4iJ  1  t 

j  =  1,2,...  N . } ,  on  S : 


where 


T.  (t) 

1  £ 

^ti  ~  ~p~  * 


P.(t) 

1  2 

— —  i 


ti+l  "  ti 


T  (t)  J 

|  i+2  i+1 


P^t)  = 


1  =  1,2 . N 


i  =  1,2, ... ,N, 


f°r  t.  <_  t  I  ti+1 


for  t.+1  <  t  <  t1+2 


for  t  elsewhere 


for  t .  1  t  <  t1+1 


for  t  elsewhere 


(3-5) 


(3-6) 


(3-7) 


(3-8) 


pi  ‘  \  (l\  +  'W 


(3-9) 


N3  -  2 


for  simply  connected 


for  doubly  connected 


(3-10) 


n3  -  1 


for  simply  connected 


for  S.j  doubly  connected 


(3-11) 


Note  that  when  is  doubly  connected,  as  shown  in  Fig.  8(a),  we  use: 


\  +  l’  ZN.+1)  ”  (pi>  z])  (pN7+2’  zN,+2^  =  (P2’  Z2)  (3-12) 

4  4  /.  /■ 


4  '  *  1  4 
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The  symmetric  product,  <a,  b>  of  two  vector  functions  is  defined  as: 


<a,  b> 


j  a  •  b  p  dt 

r 


(3-13) 


3.2.  Matrix  Equation  for  the  Nonmodal  Formulation 

The  moment  method  [17]  is  used  to  reduce  the  set  of  linear 
operator  equations  (2-76)  to  (2-78)  to  a  matrix  equation.  To  do  this, 
we  first  expand  our  unknowns  M^,  and  J11  as  linear  combinations  of 
the  basis  functions  defined  in  the  last  section: 


'"it  l<t> 

=  y  vn  u  +  y  vn  u 

-1  ^  Itq  -t,q+N2-l  ^  l<t>q  -<J),q+N2-l 


(3-14) 


m2t  W2<p 

m"  =  y  v"  u  +  l  v”  u. 

2  q=l  2tq  ~~tq  q=l  2(*>q  ^q 


(3-15) 


y=  I  in  u  +  y  i"  u, 

b-  Jl-.  tq  -tq  t.  (J>q  ^q 


(3-16) 


q=l 


q=l 


where 


"it  =  N3  N2_1  ’  m14>  N3  N2  ’  m2t  =  Ni-2,  m2iJ)  “  Nl~1  (3~17^ 

Note  that  the  boundary  conditions  that  the  components  of  m"  and  m”  normal 
to  the  ed?,e  of  the  apertures  be  zero  are  satisfied  by  (3-14)  and  (3-15). 
Equations  C'-14)  to  (3-16)  are  substituted  into  (2-78).  Next,  we  test 
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(2-76)  with  every  element  in  the  set  {u^,  u^  |  i  ■  ^,...,^-2; 
j  =  N2,  •••»  Nj-l).  In  other  words,  we  require  that  the  symmetric 
product  of  both  sides  of  (2-76)  with  every  element  in  the  above  set  be 
equal.  Similarly,  (2-77)  is  tested  with  the  set  {ij^,  u^j | i®l,2, . . . ,N^-2; 

1=1,2 . N< 


j  ■ 

1,2,... 

,  N1-l) 

and  (2- 

•78)  is  tested  with 

the  set 

{Uti> 

j  = 

1,2,... 

*  V* 

This  moment  method 

procedure  results  in 

tion 

'zan 

tt 

7an 

tip 

0 

0 

u1 

tt 

u1  1 

%  j 

V 

zan 

4>t 

zan 

H 

0 

0 

u1 

U<M> 

V14> 

0 

0 

£ 

Zc" 

1 4> 

u2 

tt 

0 

0 

cn 

<M> 

u*t 

ylbn 

tt 

lbn 

t<f> 

2bn 

tt 

2bn 

t<J> 

£ 

-bn 

t<(> 

n 

lbn 

.  V 

lbn 

2bn 

V 

y2bn 

zbn 

. 

i" 

L  *  J 

where  the  Y 

's,  Z's 

and  U ' s 

1  are  submatrices. 

Their 

elements 

row 

and  jth 

column 

are  defined  as: 

O 

pq 

ij 

%,i+N2-l' 

Len(u 
a  — q, 

j+N2- 

1>> 

p»q 

=  t ,  tj> 

i  = 

1,2,. 

’  mlp 

3  = 

1,2,. 

•  •  - 

'  miq 

O 

pq 

ij 

<u  Len(u  .)> 

-pi’  c  -qj 

-n 

Pt 


-n 


0 

0 

0 

0 


(3-18) 


(3-19) 


p,  q  -  t,<J> 

i  -  1,  2, . . . ,  m 


2p 


j  -  1,2, ,  m 


2q 


(3-20) 
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2  'c 

[U  ]  ,  =  -z —  <u  , ,  n  *  u  ,> 

pq  ij  2nb  -pi’  -  — qj 


— 1 1  (6ii + 

4nb  Pj 


if  p-t,  q-4> 


-n  A 

— £-r~  +  6,  i  ,  +  <$  <5,  N  )  if  P“4>  .  q*t 

4n.pt  1J  i_1’j  11  j’N4 

b  l 


L  o 


if  p  =  q 


p,q  =  t,<f 

i  =  1,2,...,  m2p 

j  =  1,2 . Nq  (3-25) 

Note  that  explicit  formulas  are  obtained  for  elements  of  the  U's  because 
the  evaluation  is  very  easy.  The  Kronecker  delta  is  used  in  (3-24)  and 
(3-25).  The  evaluation  of  (3-19)  to  (3-23)  requires  complicated  double 
integration.  Numerical  methods  a^e  used  in  the  computation.  The  evalu¬ 
ation  of  the  g's  and  G's  in  (2-69)  to  (2-72)  is  important  because  of  the 
special  care  required  in  handling  the  singularities  of  the  integrands. 

The  singular  part,  or  a  function  that  has  the  same  singularity,  is  inte¬ 
grated  analytically  for  each  integrand  that  has  a  singularity.  The  dif¬ 
ference,  which  is  regular,  is  then  integrated  numerically.  The  details 
are  explained  by  Mautz  and  Harrington  in  [18],  [19],  and  [20]. 

Vn  ,  v"  v"  ,  V"  ,  Tn  and  7^  are  column  matrices  defined  as 

It’  14>  2t  24>  t  <)> 
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V? 

is 


isl 


is2 


ism 


is 


(3-26) 


i  =  1,2 
s  =  t,  <j> 


in  = 

s 


si 

[\ 

s2 


sN 


t,4> 


(3-27) 


Their  elements  are  the  unknown  coefficients  of  the  expansions  in  (3-14) 

to  (3-16).  pn  and  p”  are  column  matrices  of  lengths  m.  and  m, ,  re- 
t  9  It  lep 

spectively.  Their  elements  are  defined  as 


Psj  ^a  — s.j+^-l’  — i 


s  =  t ,  <f) 

j  =  1,2, ... ,m 


Is 


A  more  detailed  examination  of  the  "excitation  matrix"  defined  in 
(3-28)  is  discussed  in  a  later  section. 


(3-28) 


3.3  Matrix  Equation  for  the  Modal  Formulation 

To  discuss  the  modal  formulation,  we  first  note  that  the  complete 
systems  of  equations  for  this  formulation  consists  of  the  operator  equa¬ 
tions  (2-79)  and  (2-80)  plus  the  equations  that  set  the  conditions  on  the 
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magnetic  currents,  when  required,  as  summarized  in  Table  1.  When  the 


,  wn 


operator  L,  is  used,  a  finite  number  of  modes  are  used  instead  of 

D 

an  infinite  number  of  them.  The  unknown  magnetic  curretns  and 
are  expanded  as  in  the  last  section: 

mlt  ml<|> 

Mn  =  y  vn  u  +  y  vn  u 

-1  i£1  vlti  ^t,i+N2-l  1f1  V10i  ^,i+N2-l 


(3-29) 


m2t  m2<j) 

‘  Jj  v5ti  a.,,  +  J1 


(3-30) 


We  substitute  (3-29)  and  (3-30)  into  (2-79)  and  (2-80).  Equations  (2-79) 
and  (2-80)  are  then  tested  by  {u^,  u^  |  i  =  N2 , .  . . ,  N^-2;  j  «=  N2,...,N.j-l} 
and  (utl,  u^j  |  i  =  1,2,...,  t^-2;  j  =  1,2,...,  1^-1} ,  respectively,  using 
the  symmetric  product  defined  in  (3-13).  The  resulting  equations  can 
be  written  in  the  following  form: 


“an 

‘5 +  C1 

bnl2 

Wtt 

bnl2 

t<t> 

rt 

1 

ztt 

zan 

V 

+  w'?”1’ 

4>t 

an  bnll 

4>ch  '<M> 

wbn!2 

4*t 

wbn12 

4>0 

x" 
i  <P 

wbn21 

tt 

wb"21 

t<J> 

7cn 

+  wbf2 

7Cn 

bn  2  2 

Vn 

Ztt 

Jtcf> 

X2t 

wbn21 

<f>t 

bn21 

zf 

4>t 

+  wbn22 
Ifrt 

zcn 

<M> 

+  wbn22 

vn 

x2ch_ 

n« 


-n 

V14> 


vn 

2t 


vn 

2<t> 


cn 


-n 

Pt 


-n 

P 


0 

0 


(3-31) 


where  the  Z's  and  V's  and  p's  are  defined  the  same  as  in  (3-26)  and 
(3-28).  The  elements  in  the  W's  are  defined  in  the  following: 
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^a  )  r  ^h  ,  n.  n  -mlD  “rnla"* 

2j7T  J  k~  C°t(knmd)  Xnm  xnm 

b  L  m=l  nm 

sink  d^O 
nm 


Ne  k' 


,  r  mn  -elp  -elq 

+  )  —  cot  (k  d)  x  rx  M 

S  k,  nm  nm  nm 

m=l  d 


sink'  d^O 
nm 


+  cot(k^d)  x^EM  x™  (if  n=0,  R^^O,  sink^d^O 


and  p=q*(J>) 


+  t-~ r  xe,lp  xe,lq  (if  there  exists  a  k*.  =  0) 
k,  d  nk  nk  nk 


P,  q=t,<|> 


(3-32) 


N  k  + 

l  —  cot(k  d)xm2pxm2q 

k  nm  nm  nm 

m-1  nm 

sink  d^O 
nm 


N 

+  I 


-rp^  cot(k'  d)  xe2p  xe2q 
m=l  k  nm  nm  nm 

sink'  d^O 
nm 


+  cot(k^d)  x^™  n=0,  R^^O,  sink^d^O 

and  p=q=4>) 


+  — r  x  p  x  .  q  (if  there  exists  a  k  =  0) 
k,d  nk  nk  nk 

D 

p,  q  =  t »<t> 


(3-33) 
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bnl2  m  _ a_ 

pq  *  2jn. 


c  _ b  ,  ,  -mlp  -m2q 

)  : -  csck  d  x  x 

,  k  run  nm  nm 

m=l  run 


sink  d^O 
nm 


Ne  k’  .  ,  + 

+  l  esc  k'  d  xelp  xe2q 

k,  run  nm  nm 

m=l  d 


sink'  #0 
nm 


+  csck,  d  x?EM  x^EM  (if  n=0,R.  #0,sink  d^O,  and  p=q=<)>) 
b  i  i  in  b 


+  t — j  x  ,  x  (if  there  exists  a  k '  =  0), 

k  d  nk  nk  nk 

b  J 

p,q  =  t  ,(J) 


(3-34) 


,,bn21  ^c  v  b  ,  -m2p  -mlq 

W  =  -r - - )  —  csck  d  x  x 

pq  2jn,  k  nm  nm  nm 

b  m-1  nm 

sink  d^O 
nm 

Ne  k'  „  ,  + 

T  ~~~  c sc  k '  d  xe2p  i6lq 
'■  k,  nm  nm  nm 

m=l  b 


sink'  d^O 
nm 

_tfm  _tfm 

+  esc  k^d  x^  (if  n=0,  R^^O,  sink^d^O  and  p=q=(j)) 


- — j  x  x  ,q  (if  there  exists  a  k'  =  0) 

k,  d  nk  nk  nk 

b 


k^d  x™j^P  x™j^q  there  exists  a  k  =  0) 


p,q  =  t ,  cf) 


(3-35) 
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,  -mlp  -elp 

where  x  ,  x  , ....  etc.,  are  column  matrices.  Their  elements  are 
nm  nm 


defined  as 


m 


, -ml t . 

Xnm  i  ^t,i+N„-l’  — nm'”  •ln  J 


put  Ti+N  -1 

I  - „»  dp 

p  Mnm 


in 


(3-36) 


i=l , 2, . . .  ,m. 


It 


r-ml4>,  .  m 

Xnm  i  “0,  i+N.-l  ’  --Tim' 


R 

■out 

p 

i+N2-l 

Pi+N2-1 

i  ,'dP 

nm 


(3-37) 


i-1 , 2 , . . . ,m 


l<t> 


f  It . 
l"'nm  ]i 


I  - 

l^t ,  i+N^-1  ’  Tim 


R 

rout 


in 


T  .  . . ,  ,  J-  dp 

l+N  ^  - 1  nm 


(3-38) 


i  =  1,2 . m 


It 


r  ~el(j) , 

x  1, 

nm  i 


■■^i+N,-!'  ^nm 


R  P  . ,  „  , 

rout  1+N„~1 
.  1  ,  6  . 

"  jn  T- .  '!>  dp 

1  p+  nm 

R.  i+N-1 

in  2 


(3-39) 


i  =  1  2 , . . . ,m 


If- 


[x 


TKM. 

1  i 


^a+N2-r  - 


TEM 


'l+N2-l 


+  — 

^  i+N-1 

*-  \ 


In  (R  / R.  ) 
out  in 


(3-40) 


i  =  1 , 2 , .  .  .  ,m 


1* 


,  -m2 1 ,  „  .  m 

x  1  .  =  -^ii  . ,  h  =  -In 
nm  i  — ti  -nm  1 


R 

rout 


i  ,  m  , 

—  >1  dp 

p  nm 


(3-41) 


in 


i  =  1,2, ... ,m 


2t 
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t\  i  rout  P,  , 

r-m2ip.  .  ,ni  1  .m'  , 

[x  ]  =  <u  ,  h  >  =  — -  ip  pd 

nm  i  —  <pl  — nm  J  +  nm 


R,  Pi 
in 


1  l,2,...,m 


2<P 


(3-42) 


,-e2t,  ,  e  _ 

[x  ]  =  <u ,  h  >  =  - 

nm  l  t  i  nm 


R 

rout 


Ti  C  dp 

i  nm 


xn 


i  =  1,2, ...  ,m 


2t 


(3-43) 


f  — e24> ,  _  .  ,  e 

Lx  .  =  *-u .  . ,  h  >  =  -  in 
nm  l  -tpl  — nm 


R 

out  P . 

i  ,  e  . 

—  *nm  dp 

R.  Pi 
in 


(3-44) 


x  1 , 2 ,  .  .  .  ,m 


2<P 


,-TEM.  ^  ,  TEM 

X2  h  =  %i*  *  "  = 


A. 

_ i _ 

3*  fin  (R  /R7T 
i  v  out  in 


(3-45) 


i  =  1,2, ...  ,m 


2 <j> 


Km’  ‘Cn’  Ti  and  Pi  are  defined  in  (2-21) ,  (2-22),  (3-7)  and  (3-8).  The 

superscript  "+"  on  a  matrix  denotes  its  complex  conjugate  transpose,  n”1 
0 

and  N  are  the  numbers  of  TM  and  TE  modes  used  in  the  approximation  of 

L  .  We  normally  choose  them  such  that  ip1"  (p)  and  (p)  have  a  number 
D  nm  nm 

of  oscillations  in  any  given  subsection. 

The  X's  and  the  column  matrix  Cn  exist  only  when  at  least  one 

"cavity  resonance"  occurs.  By  "cavity  resonance,"  we  mean  any  of  the 

situations  listed  in  Tabic  1.  The  elements  of  the  column  matrix  Cn  are 

the  unknown  coefficients  C  ,  B  ,  C  ,  and  D  when  their  corresponding 

o  nm  nm  nm 

resonance  occurs.  If  such  an  unknown  coefficient  exists  as  the  ith  ele- 


ment  of  C  ,  then  there  exists  a  corresponding  ith  column  in  each  of 


the  X^  and  .  This  corresponding  column  is  of  the 

following  form: 


’  n 

a 

-it 

X 

Y 

na 

-14> 

X 

2j 

n 

C 

cos 

kd 

-2t 

X 

n 

cos 

kd 

-20 

c 


where  the  x's  and  Y  and  k  used  in  (3-46)  depend  on  the  particular 
"cavity  resonance"  under  consideration.  The  correspondence  between 
the  unknown  coefficients  and  the  x's  and  Y  and  k  in  (3-46)  are  sum¬ 
marized  in  Table  2. 


Coefficient 

Y 

k 

-is  i=1’2 
s=t,4> 

-TEM  for  s=<j) 

C 

i/nb 

K 

Xi 

o 

0 

0  for  s=t 

B 

l  /  n. 

k 

-m 

X 

nm 

b 

nm 

nm 

-m 

c 

Y 

k 

X 

nm 

nm 

nm 

nm 

-e 

D 

Y' 

k' 

X 

nm 

nm 

nm 

nm 

Table  2.  Correspondence  between  x's  and  Y  and  k  in  (j-46)  and  the 
coefficients  associated  with  cavity  resonances. 

Next  we  consider  the  conditions  on  the  magnetic  currents  when 
cavity  resonances  occur.  Substituting  (3-29)  and  (3-30)  into  the 
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conditions  listed  in  Table  1,  we  find 


n  +  n  + 

—  x?  —  xa 

na  it  na  1* 


2t 


2<t> 


1<P 

vn 

2t 


Vn 

LV2<p 


(3-47) 


Combining  (3-38)  and  (3-54),  we  obtain,  for  the  modal  formulation. 


zan  +  wbnll 
tt  tt 

,an  bnll 

<f>t  4>t 


w 


bn21 

tt 

bn  21 


za"  +  wb"n 

t(J)  t<J) 
2“  +  «bn11 


”+t 

+ 
„n 
V  X, 
ca  It 


W 


w 


H 


w 


w 


bn21 

t4> 

,bn21 


bnl2 

tt 

bnl2 

<J>t 

bn22 


zcn  +  W 

tt  tt 

z-  +  wbn22 

<pt  <pt 


v  x" 

ca  10 


xn 

*2 1 


W 


bnl2 

t<J> 

bnl2 


W 


*5 +  2 

+  “H22 


2<D 


slt 

1<P 

‘2* 


’It 

Vn 

14> 

vn 

2t 

Vn 

20 


cn 


-n 

pt 


-n 

p* 


0 

0 

0 

(3-48) 


where 


ca 


n  /n 

c  a 


(3-49) 


3.4.  Far  Field  Measurement  and  Plane  Wave  Excitation 


The  modal  and  non modal  formulation  discussed  in  previous  sections 
share  the  same  excitation  matrix  defined  in  (3-28).  The  values  of  the 
matrix  elements  depend  upon  the  specific  type  of  incident  field  under  con¬ 
sideration.  In  this  section  we  consider  plane  wave  incidence.  This 
plane  wave  can  be  considered  as  the  field  radiated  by  an  electric  dipole 
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located  so  far  away  from  the  aperture  region  such  that  the  distance 
between  the  dipole  and  the  aperture  region  is  much  greater  than  the 
linear  dimension  of  the  aperture.  Therefore,  from  (3-28),  the  compu¬ 
tation  of  the  excitation  matrix  involves  calculating  the  reaction 
[16,  Sec.  7-7]  between  the  magnetic  field  due  to  an  electric  dipole  and 
the  basis  functions  of  the  magnetic  current  in  the  aperture  S^.  Further¬ 
more,  from  reciprocity  [16,  Sec.  3-8]  the  problem  can  be  thought  of  as 
calculating  the  reaction  between  the  electric  field  due  to  the  basis 
functions  of  the  magnetic  current  and  the  distant  electric  dipole.  This 
situation  is  very  similar  to  the  one  we  have  when  the  far  field  on  the 
transmitted  side  is  considered.  In  that  situation,  each  of  the  two  com¬ 
ponents  (tangential  to  the  radiation  sphere)  of  the  radiation  electric 
field  can  be  thought  of  as  the  reaction  between  the  electric  field  due 
to  a  linear  combination  of  the  basis  functions  of  the  magnetic  current 
and  a  distant  electric  dipole  of  unit  magnitude,  pointing  in  the 
appropriate  direction.  Therefore,  the  far  field  measurement  and  the 
plane  wave  excitation  share  basically  the  same  analysis. 

First,  let  us  consider  the  plane  wave  excitation.  The  plane 
wave  considered  here  is  of  either  the  O-polarization  or  the  <}>-polari- 
zation.  The  electric  and  magnetic  fields  of  these  two  types  are  of 
the  following  form: 


— i  o  -i 


H°  -  --  K. 
-r  r,a  -i 


O-polarization 


(3-50) 


-i 
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where 


-jk,  •  r 


—  1  ■■  C  x 

n  -i  h 

a 


<j>-polarization 


h  =  kaki  , 


(3-51) 


(3-52) 


0^  and  £  are  the  unit  vectors  in  the  spherical  coordinate  system 
at  the  point,  (r^,  0^,  <J>^),  where  the  distant  incident  source  is 
located.  This  situation  is  shown  in  Fig.  9.  A  general  plane  wave 
is  a  linear  combination  of  (3-50)  and  (3-51).  Substituting  (2-81), 
(3-50)  and  (3-51)  into  (3-28),  we  obtain 


nq  o 
3sj  ~  2 it 


Eq  rN 


3 


2tt 

f 


dtP  j 


0 


d<t»(s  •  v  £1)e 

S  =  t  ,  0 

q  =  0 » <J> 

j  =  l,2,...,mls 


“jkj  *  r-jn<t> 


(3-53) 


The  extra  superscript  q  denotes  the  type  of  polarization  under  con¬ 
sideration.  Substituting  (3-5),  (3-6)  and  (3-52)  into  (3-53),  we 
obtain 


-no 


o  ~inh 

E  e 


2"k 


n 


_  0  -i^i 

-nO  _  Eo  £ _  -<H 

^(f)  2?Tk  n 


E*  c“jn$i 
-n<t>  _  22 

Pt  "  2uk  Rn 
a 

-F*  e  1 

-n<t>  =  _o _ -<J)0 

2iik  n 


(3-54) 
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where  the  elements  of  the  column  matrices,  R's,  are  defined  as 


T1t0  ,  n+1  a 

R  .  •  j  irk  cos  8.  , 
nj  a  i  J 


Tj+N2-l(t)[Jn+l(kaP  Sin  V 


-  J  ,  (k.  p  sin  0  )  ]  sin  V  dt 
n-1  a  i 


„<|>0  .n  .  . 

R  ,  =  -  j  tt  k  cos  0. 
nj  a  i 


C-N3  Pj+N2-l(t) 

%  ‘W1 


[Jn+l(kaP  Sln  9i) 


+  J  (k  p  sin  0 . ) ]  pdt 
n-I  a  l 


„td>  ,n  , 

R  ,=  ]  TTk 
nj  a 


Tj+N2-l(t)  [Jn+l(kaP  Sin  0i> 


+  J  (k  p  sin  0.)]s invdt 
n-1  a  i 


_<b(t)  .  n+1 

RT .  =  1  tt  k 
nj  a 


t„  P..  ,(t) 

-N3  j+n2-lv 


cn2 


[Jn+l(kaP  Sin  8i) 


-  J  ,  (k  p  sin  0  . )  ]  pdt 
n-1  a  l 


(3-55) 


The  J  's  denote  the  Bessel  functions  of  the  first  kind.  The  definition 
n 

of  the  R's  is  to  link  quantities  in  this  work  to  those  in  [21]. 

Next,  we  consider  the  far  field  pattern  on  the  transmitted  side. 

From  the  equivalent  situation  discussed  in  section  2.2.,  we  know  that  the 

far  field  is  due  to  2M„  in  S„,  radiated  into  the  unbounded  medium  (e  ,  y  ) . 

i  c  c 


By  reciprocity,  we  have 
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E  •  U  =  L“(U)  •  2M.(r’)da' 

— c  —  c  —  — 2  — 


(3-56) 


where  l£  is  an  electric  dipole  situated  at  the  point  (r,0,<j>)  distant 
from  the  aperture  in  region  C.  £  is  either  _§  or  £  since  it  is  along 
the  directions  of  these  two  unit  vectors  that  the  far  field  is  to  be 
computed.  This  situation  is  shown  in  Fig.  10.  Notice  that  we  have 
moved  the  origin  to  the  center  of  S2  for  convenience.  Since  either 
0^  or  ^  is  tangential  to  the  radiation  sphere  of  12.  in  the  vicinity  of 


S^.  we  can  write 


-jk  r 

,  jk  (r  x  I?)  e  c 

Lh(lJ)  -  c  ~  ~ - 

c  —  4Tir 


£  =  0  ,c 


(3-57) 


From  (2-62b)  and  (2-62c) ,  we  obtain 


M  =  }  [Mn  £  +  MnJ>]ejn(f> 

~P  n=_,„  Pt_  P't’* 


P  =  1,2 


Substituting  (3-57),  (3-58)  and  (3-15)  into  (3-56),  we  obtain 


E  •  £  =  E 
— cr  —  cr 


(3-58) 


e"jkcr  m  Nx-2  tN  2ir 

- l[  l  V»  1  dt’T  (t ’ )  d^’i*  •(£ x  £)ej-*- +Jn4> 

‘^Trr  n=-i,>  i  =  i  ^*-3  J 

J  tx  0 


“r1  [N  p.(t')  ?’  .  ,  7 

+  ):  v“.  1  dfp’  -l  .{«!).-«•  £+j°*J 

j'1  tj  °j  0 


(3-59) 


Here  the  extra  subscript  r  denotes  radiation  field  and  the  superscript 
£  denotes  the  component  (0  or  4>)  of  the  field.  The  integrals  in  (3-59) 
are  very  similar  to  those  in  (3-53).  It  is  not  difficult  to  find  that 
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E9 

cr 

-jkcr 

-jk  e  00  .  . 

c  y  jn<|> 

27rr  1 

n=-oo 

’  -Rt4> 
n 

n 

vn  ' 

2t 

E* 

cr_ 

n 

g4>0 

n 

vn 

[>_ 

where  Che  g's  are  Che  cranspose  of  Che  R's  defined  as  in  (3-55) 
wich  (kg,  61,  tN  ,  tN  )  changed  Co  (kc>  0,  t^  CN  )  and  the  sub¬ 
script  j+N^-l  on  T,  P  and  p+  changed  to  j.  The  lengths  of  the 
matrices  also  change  from  m^,  m^  to  and  accordingly.  Note 
sin  v  is  -1  in  and  is  1  in  S2. 

It  is  appropriate  to  discuss  the  symmetry  with  respect  to  n 
at  this  time.  From  (3-56)  it  is  apparent  that 


*  -t0 

R. 

^ ' 

r-Rte 

-R^l 

-n 

-n 

n 

n 

-Rte 

R^ 

-R*9 

R* 

L  -n 

-n  „ 

_  n 

n 

From  (3-19)  to  (3-23)  and  (2-67)  to  (2-73),  we  find 


■  nt,-n 
^tt 

t<t> 

-  7<tn 
tt 

an- 
_  tCf) 

f»,-n 

Lit 

7™,-n 

H  J 

an 

L“V 

7an 

a  =  a,b,c 


(3-61) 


(3-62) 


and 
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Yjb,-n 

tt 

Yjb,-n 

t<p 

II 

_Yjbn 

tt 

Yjbn 

t<j> 

Y  j  b ,  -n 
_  <pt 

Y j  b , -n 

<M> 

Y  jbn 

L  V 

Yjbn 

<M> 

j  -  1,2 


(3-63) 


The  elements  of  the  U's  do  not  depend  on  n,  as  shown  in  (3-24)  and 
(3-25).  Therefore,  from  (3-54)  and  the  above  discussion,  we  find 


. 

V_n<^  = 

jt 

Yn(t> 

jt 

j  -1,2 

II 

CD 

C  u 

1  ‘n 
l> 

^  ■ 

v"n^  = 

j<}> 

_  vn0 
jt 

j  =  1,2 

(3-64) 
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Y~ n4»  _ 
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r  ■ 

-nfi 

> 

7“n<t>  _ 

Yn<f> 

4> 

Here  the  extra  superscripts  <3  and  <j>  denote  the  solution  of  (3-18)  for 
the  particular  polarization  of  the  incidence.  From  (3-14),  (3-15), 
(3-58)  and  (3-64) ,  we  obtain 


It 


t=t .  ,  M 
i+N„ 


nO 

<v,  y 

v  Iti 

=  2J  l  n  '  '  sin  n 
n=l  1  i+N„ 


i-1.2,, 


,o0 


ot,  V 


n0 


M 


+ 


=  +2  £  cos  n4>  i*l,2,, 


t_ti+N9-J  Pi+N2-1  0=1  Pi+N2-1 


2t 


,,n0 

on  y 

v  2ti  . 

=  21  >  - -  sin  n 


i=l,2. 


’“it 

,ml4) 

,m2t 


oo 


:i+N2  "14*2  n=1  "14M2- 


OO  V  ^ 

+  =  ^  l  5111  n<)) 

t-ti+N2-l  n  1  Pi+N2-1 

vo4>.  -  vn*. 

=  — -  +  2  l  cos  n <fi 

2t  t-t1+1  Pi+1  n-1  Pi.+1 


i  in . 


i  i  >  2  j . . .  ( m 


+ =  2j  l 


sin  n0 


n=l  p. 


i-1 . 2 f . » * ^2^  ( 


The  extra  superscript  on  M's  again  denotes  the  polarization  of  the 


incident  wave.  The  electric  current  J  exhibits  similar  behavior. 

0  0 

J  and  J,  have  the  same  angular  behavior  in  each  mode  as  that  of 
t  <p 

M?  and  M?..  and  have  the  same  angular  behavior  in  each  mode  as 
It  l<p  t  (p 

that  of  M?.  and  MP  : 

It  lcp 


06  Tn0 

I  .  m  I  . 

ti  .  „  r  ti 

=  -  +2  > - cos  n® 

D 

,t=t  Mi+1  n=l  pi+] 


i  =  1,2,...,  N 


—  »»  W 

CO  J 

V*  +  =  2j  J  -fi  Sin  n<^ 
u  v  t=ti  n=l 


i  =  1,2, ...,N, 


Vtjir  Mai 
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n 


fH 

k  • 


d> 

n  r 

Vt 


n  j+ 

Vet. 


f»  I 


=  2j  l 


ncf) 

ti 


t=t 


i+1 


n=l  Pi+1 


sin  n0 


_o0 


rn(j> 


00  I 

=  — +2  i  cos  n<j) 

n=l 


i  =  1,2,. ...N 


i  =  1,2,. ..,N 


t"ti  Pi 


0 

(3-66) 


From  (3-60)  and  (3-64),  we  obtain 
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The  first  superscript  of  E  denotes  the  component  of  the  far  field 
and  the  second  superscript  denotes  the  polarization  of  the  incidence. 
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Chapter  4 

EQUIVALENT  CIRCUIT  AND  LOW  FREQUENCY  APPROXIMATION 
OF  THE  NARROW  ANNULAR  SLOT 

In  this  chapter,  the  narrow  annular  slot  shown  in  Fig.  11  is 

considered.  Now  region  b  inside  the  thick  screen  is  a  coaxial  region. 

For  simplicity,  regions  a  and  c  are  assumed  to  be  air  filled  with 

k  =  k  =  k  and  n  =  0  =  n  •  The  word  "narrow"  means  the  following: 

a  c  o  a  c  o 


Conductor 


Fig.  11.  A  narrow  annular  slot  in  a  thick  screen,  cross  section  in 
the  y-z  plane. 
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imr 


nm 


om 


(1-Y) 

mlT 

(1-Y) 

mu 


n,ra+l  (1-Y) 


n  *  0,1,2, . . 

m  =  1,2, .. . 

in  1,2, ... , 

n  =  1,2, .. . 
in  1,2,... 


nl 


n  —  1,2, ..< 


(4-5) 


It  is  then  straightforward  to  see  that  if  (4-1)  and  (4-2)  are  true, 
and  if 

Vout  <  xn  =  1  •  (4'6 

then  the  TEM  mode  is  the  only  mode  that  propagates.  Therefore,  the 
power  associated  with  the  TEM  mode  is  of  special  interest  since,  when 
(4-6)  is  true,  it  can  be  responsible  for  most  of  the  power  transmitted. 

In  the  following  sections,  we  develop  an  equivalent  circuit 
for  the  problem  based  on  the  analysis  of  the  TEM  mode.  Resonant  be¬ 
havior  of  the  power  transmission  is  observed  and  the  electric  polari¬ 
zability  is  discussed  for  a  small  and  narrow  annular  slot  in  the  thick 
screen. 


4.1.  Equivalent  Circuit 

The  time  averaged  power  transmitted  through  the  aperture,  de¬ 
noted  Ptrans»  can  be  obtained  by  integrating  the  complex  Poynting 
ector  over  the  aperture  surface  S2-  Therefore,  we  have 


trans 


=  Re 


H.  •  z  da 
- b  ~ 


=  -  Re 


M,  •  H,  da 
—2  — b 


(4-7) 
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Substituting  (2-43)  into  4-7)  and  retaining  only  the  contribution  from 
the  TEM  mode,  we  obtain 


P 


o 

trans 


-  Re(V2  I*) 


(4-8) 


where  the  superscript  "o"  denotes  the  contribution  from  the  TEM  mode, 
and  the  quantities  on  the  right  hand  side  of  (4-8)  are  defined  in  the 
following : 


V 

1 


-icsc  (kbd)VL 


fh\  (h 


TEM 


(kbd)v2 


R 

■out 

M°  dp 
J  <p 

R. 

in 


j  =  1,2 


(4-9) 


(4-10) 


Equations  (4-9)  and  (4-10)  define  three  of  the  four  parameters  of  a 
two-port  network  that  we  introduce  later.  The  fourth  is  defined  as 


I ,  =  ~  J  cot  (k  d)V  -  esc  (k  d)V, 
in,  b  i  n.  b  i 


(4-11) 


We  now  proceed  to  develop  the  equivalent  circuit  of  the  problem, 
starting  with  the  operator  equations  (2-79)  and  (2-80).  Substituting 
(2-65),  (2-67a)  to  (2-67d)  and  (2-75)  into  (2-79)  and  2-80)  for  n  =  0, 
and  separating  the  the  4>  and  t  components,  we  obtain,  for  the  ^-component , 


Jko& 


out 


-jk  R 

o 


dp'n'M  (p')  ~°S-  - - —  d0  -  — J -  [cot(k  d)V  +csc(k  d)V  ]h 

1(t>  R  /2tt-  bl  b  2  — 


TEM 


R, 

in 


j  l  Y  [cot(k  d)(h™  ,M°)  +  csc(k  d)(hm  ,M°)]h' 
■,  om  om  om  —1  om  — om  — 2 

m=  i 


= 


m  .,o,  ,  ^m 
— om 


(4-12) 


2  R  it  -jk  R 

e  g)  rout  ,  o 

f-  dp’p'M^p*)  -  dcj> - J - [csc(kbd)V1  +  cot(kbd)V2]hT 

0  i  0  ^  \ 

in 


-j  £  Y  [esc (k  d) (h  ,  M1 )  +  cot(k  d)(h  ,  M.)lh 
S  om  om  -ora  —1  om  -mm  —1  — c 

m=l 


(4-13) 


where 


12  2 

R  =  /  p  +  p'  -  2pp'cos 


(4-14) 


Note  that  for  n  =  0,  the  TM  modes  have  only  the  4>-component  and  the 

TE  modes  have  only  the  p-component.  Next,  we  equate  the  inner  products 
tfm 

of  h  with  each  side  of  (4-12)  and  (4-13).  The  result  is 


AT  2  R  R 

v2  jk  rout  rout 


n„  k.  R , 

m  in 


gi  M14,(p,)p,dp' 


.  7  rjZ  f OUt 

rJ-  [cot  (k  d)V  +  esc  (k  d)V  ]  =  — — —  !  H°  (p)dp 

'b  b  1  b  2  vr-Tn>  l  ^ 

K  , 
in 


Ar  2  R  R 

v2  jk  rout  rout 

7=^-  dp  *i  M2.p(*',)p,d1' 

no  R.  R. 

in  in 


-  -i-  [esc  (k,d)V,  +  cot  (k  d)V  |  =  0 

n.  b  1  b  / 


(4-15) 


(4-16) 


where 


-jk  R 

d4) 

k  R  v 

o 


(4-17) 


'*  T'  ~rr~r  '  •  -  - 

,..Y  Am  .  ...  •**  .  ^  .^W 
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and  is  the  same  as  defined  in  (2-69).  Under  the  condition  (4-2), 
can  be  written  as 


if 

*1  k  R  1 
o  a  L 


[In 


0  2k  R 

e2 1 p  -  p1 |  !  f°  a 

8R  2 


E2(x)dx] 


2k  R 
f°  a 

+  "V  J2(x)dx 


(4-18) 


where  and  E7  are  the  Bessel  and  Weber  functions  of  the  second  order. 

The  proof  of  (4-18)  is  given  in  Appendix  A.  Substituting  (4-18)  into 

(4-15)  and  (4-16)  and  using  the  approximation  p'  =  R  ,  we  obtain 

a 


k  v 
o 


o  v 


2k  R 
o  a 


J2(x)dx  +  ^ 


2k  R 
r  o  a 


E2(x)dx]V  -  ■*'—  , 

nw/-lnY  i 


R  R  _ 

out  rout  l i 


dp 


R.  R . 
m  in 


ln(e . le--^-l-)M° 

K  8R  1<| 

a 


R 


2/2tt 


-  -f-  [cot (k  d) V  +  csc(k  d)V  ]  =  — 
nb  b  1  b  2  /-InY 


rout 


VP)dP 


(4-19) 


R. 

m 


k  w  r 
o 


2k  R 


o  a 


•J2  (x)dx  +  ^ 


2k  R 
r c>  a 


E2(x)dx]V2 - i-~L_ 

nwAlnY 


R  R  . 

out  rout  1 1  ,  | 

dp| 

R.  R.  3 

in  in 


[esc  (k  d)V  +  cot  k,  d  V  =  0 
H,  b  1  b  2 

b 


(4-20) 


We  now  examine  the  double  integral  in  (4-19) .  Changing  the  order  of 


(P')dp' 


(P')d  ’ 


integration,  we  obtain 
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where 


out 


R 

/•out 


dp 


\£=£L 


8R 


■)dp'  = 


out 


M^(p')f(p')dp'  (4-21) 


in 


in 


in 


f(p’) 


R 

rout 


ln<  8#^^ 


in 


=  w  In  — - w  +  (R  -  p  ) in  (R  -  p  )  +  (p  -  R  )ln  (p1  -  R  ) 

8R  out  out  in  in 

a 


(4-22) 


f(p')  is  an  even  function  about  p'  =  R  and  has  its  minimum  value 


2 

f  .  =  (w  In  r—  -  w)  occurring  at  p'  =  R,  and  p' 

min  HR  in 

a 

o 


R  .  Its  maximum 
out 


value  f  =  (w  In 

max  16R 


w)  occurs  at  p'  =  R  .  When  w  <<  R  ,  f(p')  is 
a  a 


almost  constant.  Therefore,  the  value  of  the  double  integral  in  (4-21) 
is  essentially  proportional  to  the  integral  of  i.e.,  V^,  ,  and  is 

insensitive  with  respect  to  the  actual  functional  form  of  M^(p'). 
Various  values  can  be  obtained  for  this  integral  by  assuming  different 
functional  forms  for  M°. (p * ) .  They  converge  to  the  same  value  as  the 


!<!> 


w 


ratio  —  becomes  very  small .  The  same  argument  is  valid  for  the  double 
K 

a 

integral  involving  in  (4-20).  Equations  (4-19)  and  (4-20)  are  now 
written  in  the  fol lowing  form: 


vh8vi  +  'Vi  +  VUV2 


=  I 


(4-23) 


Y^SV  +  Y*5  v  +  y  k  v 
2  21  1  22V2 


=  0 


(4-24) 


where 


cot  k,  d 

\  b 


esc  k.d 
nb  b 


— cot  k.  d 

n,  b 


(4-25) 


YhS  =  G  +  jB 


(4-26) 


.  2k  R 
k  w  fo  a 

G  =  2^~  Vx)dx 

°  0 


(4-27) 


rout 

!  M  .  , 
J  trial 


(p')f(p’)dp’ 


kwr,2{oRa  R , 

B  =  f~  I  E2 (x)dx  -  ^ - 

o  i  J  rout 

Mtrial(p’)dp’ 


(4-28) 


/-lny 


H^(p)dp 


(4-29) 


In  (4-28),  represents  any  particular  functional  form  that  we  choose 

for  M, ,  and  M°  .  If  we  choose  to  use  a  constant,  we  obtain 
l<p  2(J> 


k  w 

B  =  -±- 

% 


Zk  r 

’  I  f° 

.  2 


E2(x)dx  -  -  In  (~gF 


(4-30) 


If  we  choose  to  use  for  M  ,  ,  the  quasistatic  solution  for  a  narrow 

trial 


annular  slot  in  a  infinitely  thin  screen  [24],  then 
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Mtrial(p> 


/(f)2  -  (P  -  Ra)2 


(4-31) 


and  the  result  Is 


,  2k  R  _ 

k  w  r  ,  fO  a  .  2 

f-  2  E2(x)dx  -  £  ln  (M->] 

o'-  i.  a  J 


(4-32) 


Equations  (4-23)  and  (4-24)  can  be  viewed  as  the  equations  for  the 
equivalent  circuit  shown  in  Fig.  12.  This  equivalent  circuit  can  be 
used  to  analyze  the  TEM  mode  as  long  as  @t-2)  holds.  The  current  source 
can  be  computed  for  plane  wave  incidence: 


H. 

— l 


-E  -jk.  •  r 

o  i  —i  — 

— —  <■».  e 
n  -"-i 

o 


(4-33) 


Note  that  (4-33)  represents  the  0-polarized  plane  wave  defined  in 
(3-50),  since  a  <})-polar  ized  plane  wave  does  not  excite  the  TEM  mode. 
Substituting  (4-33)  and  (2-81)  into  (4-29)  and  using  the  condition 
(4-1),  we  obtain 

-  j  2/2?  w  e'1  J,  (k  R  sin  0,) 

Ts  =  _ _ _ i_  ( 

n  /-lny 

o 

In  the  above  analysis  we  have  assumed  sin  k,d  i  0.  This  is  evident 

D 

because,  when  deriving  (4-8),  (4-12),  and  (4-13),  we  used  (2-43)  and 
(2-75)  in  the  case  where  sin  k^d  #  0.  It  is  not  difficult  to  show 
that  when  sin  k^d  =  0,  we  have  the  following  equations: 


*  * 
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hs  r 

Y  V.  +1,-1 
11s 


YhSV2  +  1X2  =  0 


(4-35) 


(4-36) 


V.  +  V.  cos  k,d  =  0 
i  l  b 


(4-37) 


r  j- 

where  and  V2  are  defined  the  same  as  in  (4-10)  while  and  I2 


are  defined  as 


i*  =-i^c 
1  \  ° 


(4-38a) 


j/2TT  cos  k  d 

ll  - - - - C 

2  \ 


(4-38b) 


C  is  the  extra  unknown  constant  associated  with  the  case  sin  k,d  =  0, 
o  b 

introduced  in  Chapter  2.  The  equivalent  circuit  for  this  case  can  be 
constructed  from  (4-35)  to  (4-37)  and  is  shown  in  Fig.  13.  The  power 
transmitted  can  be  shown  to  be 


-  Re(V2  ij  ) 


(4-39) 


for  this  case. 


4.2.  Power  Transmission  and  Resonant  Behavior 

The  equivalent  circuit  shown  in  Fig.  12  is  very  similar  to 
that  of  a  narrow  TE  slot  in  a  thick  screen  [25].  When  the  thickness 
of  the  screen  is  close  to  a  multiple  of  a  half  wavelength,  the  waveguide 
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region  can  "tune"  the  aperture  to  transfer  peak  power.  We  show  this 
by  considering  our  equivalent  circuit  in  Fig.  12.  From  the  equivalent 
circuit,  we  obtain  the  power  transmitted  as 


P°  =  -  V  I*  =  — 

trans  2L2  y. 


(4-40) 


Here  y ^  is  the  transfer  admittance  defined  as 


]  2  V„ 


(4-41) 


and  can  be  solved  from  the  circuit  as 


y12  =  -  2Yhscos  k^d  -  ^  [1  +  (nbYhs)2]sin  k^d 


(4-42) 


Equations  (4-40)  and  (4-42)  are  valid  even  when  sin  k,d  is  zero.  This 

b 

can  be  easily  shown  by  solving  the  equivalent  circuit  in  Fig.  13. 
Substituting  (4-26),  (4-27),  (4-32),  (4-34),  and  (4-42)  into  (4-40), 
we  obtain 


BttwR  J,(k  R  sin  0,)G' 
a  i  o  a  i 


trans  no  [ 2G'coskbd  -  2G' B' sink^d) ] 2 +  [2B' cosk^d +  (1+G' 2-B' 2)sinkbd]2 


where 


(4-43) 


.  2k  R 
k  w  ro  a 

G'  =  H  G  =  ■  J,(x)dx 

O  l  L 


(4-44) 


B'  =  n  B 

O 


.  r  2k  R 
k  w  f  co  a 

-H  I 


F2(x)dx  +  f  [In  -  2] 


(4-45) 
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For  a  finite  k  R  ,  under  the  condition  (4-1),  we  see  from  (4-44) 
o  a 

and  (4-45)  that 

G'  «  1  ,  B'  <<  1  (4-46) 

With  (4-46)  in  mind,  we  see  from  (4-43)  that,  when  sin  k^d  is  not  small, 

P^rans  aPProximately  proportional  to  csc^k^d  and  has  its  minimum  in 

the  neighborhood  of  d  =  (x  +  t)  •  Furthermore,  P°  is  large  in 

24b  trans 

the  neighborhood  of  d  =  ^-A,  .  Looking  at  the  denominator  of  P°  in 

2b  trans 

(4-23),  we  see  that  when  k.d  =  nir ,  the  first  term  of  the  two  is  very 

b 

insensitive  to  a  small  change  in  k^d  and  has  the  approximate  value  of 
2 

4G'  .  The  second  term  however,  is  very  sensitive  to  a  change  in  k^d, 

and  can  be  made  zero  by  precisely  choosing  our  k^d.  Therefore,  the 

maximum  of  P°  occurs  when 
trans 

tan  K  d  =  -  — - ~ - r—  (4-47) 

b  (1  +  G'  -  B'  ) 

Equation  (4-47)  is  our  condition  for  "slot  resonance"  described  above. 
Since  B'  is  small,  a  first  approximation  of  the  solution  of  (4-47)  is 


(!  -  r>Ab 


(4-48) 


The  peak  value  of  P°  can  be  found  by  substituting  (4-47)  into 
trans 

(4-43).  The  result,  after  proper  approximation,  is 


po , res 
trans 


|  E~  (  ^  2ti  w  R  (k  R  sin  0.) 
O'  a  1  o  a  i 


(4-49) 


The  superscript  "res"  that  appears  in  (4-48)  and  (4-49)  denotes  the 
slot  resonance,  and  is  to  be  distinguished  from  the  superscript  "r" 
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in  (4-35),  (4-36),  (4-38),  and  (4-39)  where  it  denotes  the  situation 
when  sin  k^d  is  exactly  zero. 

Another  quantity  of  interest  is  the  width  of  the  power  trans¬ 
mission  peak  just  discussed.  We  define  this  peak  width  as  the  dis- 
tance  between  the  two  neighboring  points  of  d  where  half  of  the 
peak  power  is  transmitted.  Therefore,  we  solve  the  equation 


P 


o 

trans 


1  p0>res 

2  trans 


(4-50) 


Substituting  (4-43)  and  (4-49)  into  (4-50)  and  using  the  approxi- 
nations 

V  =  vres  +  * 

cos  kjd  =  (-l)n  (4-51) 

sin  k,d  =  (~l)n(6  -  2B')  , 
b 

we  obtain 

6  =  +  2G’  (4-52) 
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where  P  is  defined  in  (4-7).  P.  is  the  power  that  would  be 

trans  in 

intercepted  by  the  aperture  S^  if  the  incident  plane  wave  were  normal. 
Therefore, 


2it  w  R  |  E°  1 2 

P  - - - I_ 

in  n 

o 


(4-55) 


where  E°  is  the  magnitude  of  the  incident  electric  field  of  the 
general  plane  wave  of  the  form 


ii-  «X  +  E^)e'Jti‘s 


(4-56) 


For  the  problem  discussed  in  this  chapter,  if  (4-6)  is  valid,  we  have 


P 


trans 


=  P 


o 

trans 


(4-57) 


Substituting  (4-43),  (4-55),  and  (4-57)  into  (4-54),  we  obtain 

|E°|2  4G’J?(k  R  sin  0.) 

.p  ^  _ 2. _ 1  on _ i _ _ 

|  E°  |  2  ^G'coskjd  -  2G,B,sinkbd]2  +  [aB’cosl^d  +  (l-HJ,2-B,2)sinkbd]2 

(4-58) 

Similarly,  if  (4-49)  is  used  instead  of  (4-43),  we  obtain,  for 
the  slot  resonance, 


Ki2  v 

! E° I 2g ' 


(4-59) 


4.3.  Small  Apertures  and  the  Electric  Polarizability 

If  in  addition  to  the  conditions  (4-1)  and  (4-2),  we  assume  that 
the  overall  size  of  the  aperture  is  small,  i.e., 


k  R  «  1 
o  a 


(4-60) 
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then  some  of  the  formulas  In  the  previous  sections  can  be  simplified. 
These  formulas  are  written  below  without  detailed  justification  be¬ 
cause  they  involve  mostly  simple  small  argument  approximations  of  re¬ 
lated  functions. 


7  (k  w) 

(k  R  )3 

(4-61) 

6  0 

0  a 
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In  (4-63)  to  (4-65),  E.  (0)  represents  the  z-component  of  the  incident 

iz 

electric  field  evaluated  at  the  aperture  S, .  Also,  when  R  is  small, 

1  a 

the  fields  in  region  c  can  be  considered  as  due  to  an  electric  dipole 
in  S2  in  the  presence  of  the  shorted  conducting  screen.  The  dipole 
can  be  computed  from  the  magnetic  current. 
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Substituting  (4-41)  and  (4-63)  into  (4-66),  we  obtain 
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U  =  ju>  2^ 


(4-70) 
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Chapter  5 
NUMERICAL  RESULTS 


The  computational  results  of  some  typical  examples  are  pre¬ 
sented  In  this  chapter.  Five  simple  configurations  are  used.  The 
cross  sections  of  the  apertures  in  the  y-z  plane  for  these  configu¬ 
rations  are  shown  in  Figs.  14(a)  to  14(e).  Two  narrow  annular  aper¬ 
tures  shown  in  Figs.  36(a)  and  36(b)  are  used  to  demonstrate  the 
resonant  behavior  described  in  Chapter  4.  Results  from  the  matrix 
solution  and  the  analytical  prediction  are  compared  and  discussed. 

For  the  current  plots,  the  magnetic  current  is  normalized  with  re¬ 
spect  to  the  amplitude  of  the  incident  E-field,  and  the  electric 
current  is  normalized  with  respect  to  the  amplitude  of  the  incident 
H-field.  The  horizontal  axis  represents  the  variable  t  along  the 
generating  curve.  Tick  marks  are  placed  to  show  how  subsections 
along  the  curve  are  arranged  in  the  matrix  solution.  The  (^-component 
of  either  the  magnetic  or  the  electric  current  is  represented  by  a 
combination  of  pulses,  consistent  with  the  true  expansion  functions 
in  the  solution.  The  t-component  is  represented  by  a  linear  func¬ 
tion,  constructed  by  connecting  points  representing  the  solution  at  each 
t  's  with  straight  line  segments.  This  is  an  approximation  to  the  true 
expansion  because  the  expansion  functions  for  the  t-component  are 
triangles  divided  by  p.  For  cases  to  which  the  modal  solution  also 
apply,  the  results  from  both  formulations  are  shown  for  comparison. 
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The  <j>  and  t-components  of  the  modal  solution  are  marked  with  tri¬ 
angles  and  squares,  respectively,  at  a  number  of  discrete  points. 

Figures  14(a)  and  14(b)  show  typical  cases  where  both  formulations 
are  used.  For  the  circular  aperture  in  Fig.  14(a),  normal  incidence 
is  considered.  Only  the  n  =  +  1  modes  are  needed  for  this  case,  and 
the  n  =  1  mode  currents  are  shown  for  an  incident  E-field  polarized 
in  the  -x  direction.  The  azimuth  angle  of  incidence  <{)^  is  assumed  to 
be  zero  without  loss  of  generality  because  of  the  symmetry.  The 
magnetic  current  on  the  transmitted  side  shown  in  Fig.  16  is  much 
smaller  in  magnitude  than  that  on  the  illuminated  side  shown  in  Fig.  15, 
and  the  attenuation  can  be  clearly  seen  from  the  electric  currents, 
which  represent  components  of  the  tangential  H-field,  in  Fig.  17.  The 
wall  thickness  in  this  case  greatly  reduces  power  transmission.  The 
aperture  with  a  filled  coaxial  region,  shown  in  Fig.  14(b),  exhibits 
different  transmission  characteristics.  An  oblique  wave  incident 
from  0^  =  150°  is  applied.  The  aperture  is  small  electrically,  and 
therefore  only  the  n  =  0,  +  1  modes  are  needed.  Figure  18  shows  the 
circulating  magnetic  currents  on  each  side  of  the  aperture  and  the 
t-directed  electric  current  on  the  walls,  for  n  =  0,  due  to  a  0-polarized 
incidence.  This  polarization  couples  to  the  TEM  mode  of  the  coaxial 
region.  The  effect  of  this  propagating  mode  is  seen  in  the  current  dis¬ 
tribution.  We  don't  see  the  great  attenuation  observed  in  the  previous 
case,  where  all  modes  are  evanescent.  The  currents  for  the  n  =  1  mode 
are  also  shown  in  Figs.  19  and  20.  The  attenuation  of  this  mode  is 
again  present,  and  therefore  the  TF.M  mode  transmits  most  of  the  energy. 


Fig.  15. 


M.  and  M, ,  of  the  example 
It  1<P 


in  Fig.  14(a) 


Fig.  16.  M2t  and  - JM2^  of  the 
example  in  Fig.  14(a). 
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Currents  for  the  ^-polarized  incidence  are  shown  in  Figs.  22  to  24. 

The  currents  for  n  =  1  mode  are  shown  and  they  are  very  similar  to 
that  for  8-polarization.  However,  these  currents  are  responsible 
for  most  of  the  energy  transmission  in  this  case,  since  power  trans¬ 
mitted  by  the  n  =  0  mode  is  negligible  now  that  it's  not  coupled  to 
the  TEM  mode.  The  currents  for  the  n  =  0  mode  in  this  case  are  not 
shown  for  this  reason.  Results  for  the  magnetic  currents  from  both 
formulations  are  compared  for  the  above  two  configurations.  The 
overall  agreement  is  very  good.  The  proper  edge  behavior  of  the 
magnetic  current  is  also  observed.  One  of  the  simplest  geometries 
for  which  the  modal  approach  is  not  applicable  is  shown  in  Fig.  14(c). 
The  n  =  1  mode  currents  are  shown  in  Figs.  25  to  27  for  normal  inci¬ 
dence.  The  magnetic  currents  again  exhibit  a  distribution  similar  to 
that  which  exists  in  a  small  circular  aperture  in  a  infinitesimally 
thin  screen.  The  fields  decrease  rapidly  from  one  side  to  the  other. 

It  is  interesting  to  observe  that  most  of  the  attenuation  occurred 
before  the  sharp  corner  on  the  inner  wall  of  the  aperture.  The  con¬ 
figuration  shown  in  14(d)  models  a  gasket  in  the  thick  screen.  Oblique 
incidence  of  135°  with  0-polarization  is  considered.  Figure  28  shows 
the  circulating  magnetic  currents  and  the  t-component  of  the  electric 
current  in  the  n  =  0  mode.  It  is  apparent  that,  although  our  waveguide 
region  is  no  longer  a  straight  coaxial  region,  the  existence  of  the 
center  conductor  still  enables  the  fields  to  propagate  through  the 
screen.  Also,  because  of  the  size  of  the  aperture,  there  is  propaga¬ 
tion  even  for  the  n  =  1  mode.  This  can  be  seen  from  Figs.  29  to  31. 
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Again,  the  ^-component  of  the  magnetic  currents  becomes  large  at 
both  edges  of  the  aperture  as  expected.  Our  last  example  for  cur¬ 
rents  is  shown  in  Fig.  14(e).  The  incident  wave  is  0-polarized  with 
a  135e  angle  of  incidence.  Currents  for  both  the  n  =  0  and  n  =  1 
modes  are  shown  in  Figs.  32  to  35.  The  system  is  below  cutoff  and 
the  screen  effectively  attenuates  th  field. 

The  transmission  coefficient  T  is  computed  for  the  two 
narrow  annular  slots  shown  in  Figs.  36(a)  and  36(b).  Figure  37 
shows  the  transmission  coefficient  for  the  aperture  in  Fig.  36(a) 
as  a  function  of  the  screen  thickness  d.  The  computed  results  does 
show  a  repeated  resonant  phenomenon  with  period  y.  Since  the  peaks  are 
in  general  very  narrow,  only  the  result  in  the  neighborhood  of  the  first 
resonance  is  shown.  The  solid  line  represents  the  result  from  (4-58) 
and  the  circles  represent  the  result  from  our  modal  matrix  solution. 

The  agreement  is  excellent  except  that  the  centers  of  the  peaks  are 
separated  by  a  few  ten  thousandths  of  a  wavelength.  This  can  be  ex¬ 
plained  as  follows.  From  (4-53),  the  width  of  our  resonant  peak  is 

2g  *  g 1 

approximately  -  A,  .  The  center  of  the  peak,  from  (4-48)  is  —  A 

7Tb  7T  D 

away  from  half  wavelength.  For  our  particular  problem,  G'  and  B'  can 
be  found,  from  (4-44)  and  (4-45),  to  be  roughly  0.00014  and  0.038. 
Therefore,  a  combined  error  of  one  percent  from  the  two  calculations 
of  B'  can  cause  the  peaks  to  separate  completely  from  each  other. 

and  V2,  the  voltages  defined  in  (4-10),  are  computed  from  the  matrix 
solution  also.  They  are  shown  in  Fig.  38.  The  solid  lines  are  the 
real  and  imaginary  parts  of  V^.  The  squares  and  triangles  are  the  real 
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and  Imaginary  parts  of  V^.  Note  that  the  peak  occurs  where  the  real 

parts  of  and  V ^  are  zero.  This  is  expected  because  resonance  occurs 

when  the  imaginary  part  of  the  transfer  admittance  y^  Is  zero  and  our 

current  source  is  imaginary.  Also  note  that,  since  d  is  close  to  j  , 

and  V2  are  almost  exactly  the  same  as  for  the  case  d  *  ~  shown  in 

Table  1.  Figures  39  and  40  show  the  transmission  coefficient  and  the 

voltages  for  the  case  in  Fig.  36(b).  Again,  the  agreement  is  excellent 

G  *  1 

and,  this  time,  since  -  — ,  a  comparable  percentage  error  in  B'  would 

D  ZU 

not  cause  the  pulses  to  look  completely  separated. 


1305887 X 
1240592X 


Fig.  39.  Transmission  coefficient  of  the  aperture  in  Fig.  36(b) 
as  a  function  of  the  screen  thickness. 
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Fig.  40.  Vi  and  V2,  voltages  of  the  equivalent  circuit  of  the 
aperture  in  Fig.  36(h),  normalized  with  respect  to 
X,  E°,  as  functions  of  the  screen  thickness. 
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Chapter  6 
CONCLUSION 


A  general  formulation  is  developed  for  solving  the  problem  of 
electromagnetic  transmission  through  a  rotationally  symmetric  aperture 
in  a  conducting  screen  of  finite  thickness.  The  solution  obtained  is 
in  the  form  of  the  Fourier  coefficients  (in  <p)  of  the  equivalent  cur¬ 
rents  on  the  boundary  of  the  aperture  region.  All  field  and  power 
characteristics  can  be  computed  from  these  coefficients.  The  number 
of  Fourier  modes  needed  depends  on  the  size  of  the  aperture  and  the 
nature  of  the  excitation.  Because  of  the  symmetry  of  the  problem,  each 
mode  can  be  solved  separately.  The  basic  formulation  resembles  that  of 
[10],  except  that  our  problem  is  a  three-dimensional  one.  Also,  for  the 
modal  formulation,  when  the  waveguide  region  is  of  resonant  size,  problems 
can  arise  in  forming  the  field  operator  in  this  region,  as  observed  in  [7]. 
This  problem  is  solved  in  Chapter  2  by  using  the  proper  constraints 
(Table  1)  on  the  magnetic  currents  and  introducing  new  unknowns.  Note 

should  be  made  that  the  case  k'  =0  needs  to  be  treated  carefully,  as 

run 

shown  in  Chapter  2.  However,  this  case  does  not  correspond  to  a  cavity 
resonance  While  the  procedures  for  treating  a  resonant  case  may  seem 
applicable  to  only  a  few  special  cases,  the  concept  can  be  extended  to  a 
more  general  case  because  the  magnetic  currents  on  the  aperture  faces 
can  always  be  separated  into  two  parts,  one  a  linear  combination  of  wave¬ 
guide  modes  and  the  other  orthogonal  to  them.  The  modal  solution  can 
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also  be  extended  to  treat  an  aperture  with  region  b  composed  of  cascaded 
waveguide  sections. 

In  the  low  frequency  discussion  of  Chapter  4,  an  equivalent  cir¬ 
cuit  was  developed  for  a  narrow  annular  slot.  An  effort  was  made  to  show 
the  validity  of  a  one  term  moment  solution  in  this  case,  as  done  in  [25]. 
Also,  the  assumption  was  made  that  only  the  TEM  mode  propagates.  This  limits 
the  overall  size  of  the  aperture  as  shown  by  (4-6).  If  the  mean  radius  of 
the  narrow  annular  slot  is  arbitrary,  then  there  may  be  a  propagating  mode 
for  each  n,  as  seen  in  (4-5).  Further  study  of  each  of  those  modes  is 
recommended  to  develop  the  corresponding  equivalent  circuits  for  such  cases. 
Our  analysis  of  the  TEM  mode  remains  valid  as  far  as  the  power  transmission 
associated  with  the  particular  mode  is  concerned. 
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Appendix 

PROOF  OF  EQUATION  (4-18) 

In  this  appendix,  we  prove  that 
IT 
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E2 (x)dxj 
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J2(x)dx  j 


when  W  «  R  .  R  Is  defined  as 
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;2  2 

p  +  p'  -  2pp’  cos  <j) 


The  subscript  of  k  in  (4-18)  is  left  out  for  convenience. 
Y  as  a  function  of  k,  and  write 
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Y(k)  =  Y (0)  + 

From  (A-l),  we  find 
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Since  Y(0)  ie  real,  we  have 


Im(Y) 


“M 


cos  <p  cos  (aR)d<j) 


(A-6) 


When  w  «  R  ,  we  can  obtain,  from  (A-2) , 
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R  -  2R  sin  f 
a  2 


(A— 7  ) 


Substituting  (A-7),  and  the  new  variable 


ij^  =  -  (tt  -  <Ji) 


(A-8) 


into  (A-6) ,  we  obtain 


k  tt/2 


Im(Y)  *  2 


da  cos  2ip  cos  (2aR  cos  \J;)diJ> 

J  J  a 
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(A-9) 


From  10.11.1.  of  [26],  we  obtain 


Im(Y) 


J2(x)dx 


(A-10) 


Now  we  consider  the  real  part  of  Y.  From  (A-3)  to  (A-5),  we  obtain 


Re(Y)  -  Y (0) 


* 
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dk  cos  4>  sin  (otR)d<J> 


(A-ll) 
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Substituting  (A-7)  and  (A-8)  into  (A-ll),  we  find 

k  tt/2 


Re(Y)  -  Y(0)  =  2  da  cos  2ip  sin  (2aR  cos  Tjr)d^  (A-12) 

a 
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From  10.11.3.  of  [26],  we  obtain 


Re(Y)  -  Y(0)  -  — 


2kR 
f  a 


E2(x)dx 


Finally,  we  consider  Y(0).  Substituting  (A-2)  and  (A-8)  into  (A-5 
we  obtain 
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where 


82  = 


— 4— --  V  =  1  (for  w  «  R  ) 
(P  +  c’)2 


From  17.3.1,  3,  26  of  [27],  we  find,  for  w  «  R  , 
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(A-13) 

), 

] 

(A-14) 

(A-15) 

(A-16) 


Combining  (A-10) ,  (A-13)  and  (A-16),  we  obtain  (A-l). 
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